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1. Introduction 



Flux vacua of type IIB string tlieory tiave been under detailed study for some time now 
(see [ll|fl,^3,|5|,|6|,|7|,^p|, p!0| , Pi|T2| and references therein) . Much of their interest derives from 
the fact that, in the presence of nontrivial RR and NS three-form fluxes Fs and Hs, the 
IIB theory on a Calabi-Yau orientifold manifests a computable potential for the complex 
structure moduli and the dilaton. This potential generically has isolated critical points 
for these moduli, and such vacua provide good starting points for the inclusion of further 
effects that can stabilize the remaining geometric moduli |lT3| , |T4| , p!5| , p!6[] . Similar mechanisms 
have been proposed in the other corners of the M-theory parameter space including the 
lid limit Ijl^, the heterotic theory [|18[, the II A theory [|19[, and supercritical string theory 
||20[] . One interesting feature of the flux potentials is that they give rise to a very large 
"discretuum" of vacua [|21|,^^. For this reason, it is useful (and probably necessary) to 
find statistical descriptions which can summarize the features of large classes of vacua at 
once. Such a statistical theory, which is useful for at least the simplest questions about the 
IIB Calabi-Yau flux vacua (questions such as "how many vacua are there?" and "where are 
they located on the moduli space?"), was proposed in [^,^. It has been further refined, 
applied and tested in [^|23,|2|,|9|,|g,gTl . 

In this paper, we continue this statistical study of IIB flux vacua. In addition to 
addressing the questions of vacuum degeneracy and distribution on the moduli space in 
some simple examples, we also initiate the study of two more detailed questions which may 
be relevant for models of low-energy particle physics. The first question is, what fraction of 
the IIB flux vacua are "supersymmetric" in the sense of i.e. have W = in the leading 
approximation'/S This question is important because in recent discussions of the typical 
scale of supersymmetry breaking in the landscape p^j2^ , |33| , p^ , it has been argued that 
the fraction of W = vacua may decide the issue of whether the typical SUSY breaking 
scale is high or low. A very small fraction of W = vacua would favor high scale SUSY 
breaking. This observation has motivated interesting new testable models of weak scale 
physics ||3^. On the other hand, as discussed in |3^], even a modest fraction of W = 



^ This should not be confused with the notion of supersymmetry generally used in [25| and 



[26], where the Kahler moduli are ignored and all of the no-scale vacua are called supersymmetric. 
The latter language is appropriate in many cases because, after proper inclusion of both Kahler 
moduli and nonperturbative effects, many no-scale vacua do indeed give rise to supersymmetric 
AdS vacua fU 
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vacua could lead (after further effects) to a greater number of models with low SUSY 
breaking scale. We will remain neutral about the question of preferred scale, and simply 
present results about the relative fraction of W = vacua in some tractable examples. 

The second question we address is, how frequently do discrete symmetries arise in 
flux vacua? Such symmetries are required in many models of weak scale physics, to, e.g., 
prevent a catastrophic rate of proton decay. While our simple models (which do not even 
include a standard model sector) are by no means realistic laboratories in which to face 
these detailed issues, they are a good place to begin understanding how frequently discrete 
symmetries may arise "accidentally" in the landscape. One can envision a systematic 
improvement of the realism of examples, by working through the same kinds of calculations 



in the pseudo-realistic models presented in e.g. [p6|,|37|,p8|,|39| 



Throughout the paper we use a combination of techniques to analyze the various 
models of interest. For some questions involving enumeration of the total number of 
vacua we use a continuous approximation to the fluxes, along the lines of the techniques 
developed in [p5| , |26[] . For overconstrained systems like the set of W = vacua, these 
techniques are not directly applicable, and the discrete nature of the fluxes must be taken 
into account. We develop several simple methods for dealing with such problems, which are 
more essentially number-theoretic in nature. Our approach to dealing with these problems 
relies on the algebraic structure of the periods of the compactiflcation manifold at the 
points in moduli space associated with flux vacua. In all cases we consider, these periods 
lie in a fleld extension of flnite degree over the rationals Q, so that we can reduce the 
problems of interest to a flnite set of algebraic equations over the integers. In most cases 
we have also used numerical methods to check the analytic predictions. We summarize 
our results in the flnal section. 

Our paper is organized as follows: we begin in §2 by presenting some material which 
will be generally useful in understanding our analysis. In §3-§5, we answer some of the 
questions of interest in a particular compactiflcation, with the compactiflcations arranged 
in order of increasing complexity. In §3, we study vacua in a rigid Calabi-Yau. In §4, we 
study vacua on the torus T^. We focus attention in particular on the subclass of vacua 
where the torus has the symmetric form (T^)"^, where each has the same modular 
parameter r. In §5, we describe our results (at and near a special locus in moduli space) 
for the four Calabi-Yau hypersurfaces with a single complex structure parameter. Finally, 
we close in §6 with a summary and a discussion of the implications of these results and 
promising directions for future investigation. 
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2. Taxonomy of IIB flux vacua 

The addition of fluxes can lift the moduh of string theory compactiflcations, leaving 
isolated vacua at various locations around what was the moduli space. We are interested 
in the "taxonomy" of these vacua - their general distribution as well as the incidence of 
vacua with certain special properties. 

In this section, we review the "imaginary self-dual" flux vacua in type IIB string 
theory that we will study in the rest of the paper. After describing the flxing of the dilaton 
and complex structure moduli, we briefly discuss the Kahler moduli and our philosophy 
towards their treatment. We then describe the modular symmetries that must be fixed in 
order to count inequivalent vacua, before reviewing the statistical methods applied to the 
problem of enumerating vacua by Douglas and collaborators, which primarily approximate 
the fluxes as continuous. We introduce some basic ideas from number theory and discuss 
how these ideas will help in our analysis of counting problems which are not amenable to a 
straightforward continuous approximation. Finally we motivate and describe the two types 
of "special" vacua we will be interested in counting: those with vanishing superpotential 
W = 0, and those preserving discrete symmetries in the low-energy theory. 

2.1. Basic equations 

As the properties of flux superpotentials on Calabi-Yau orientifolds in type IIB string 



theory have been reviewed many times, we will be brief. Our conventions are those of ||27|| . 

We consider a Calabi-Yau threefold M with /i2,i complex structure deformations, and 
choose a symplectic basis {A", Bi} for the 63 = 2/i2,i +2 three-cycles, a, 6 = 1, ... , /i2,i + 1, 
with dual cohomology elements such that: 



ab = S^, / /3" = -5,", / aaAf3'' = Sl (2.1) 

Fixing a normalization for the unique holomorphic three-form O, we assemble the periods 
z°' = O, = Jg^ O into a 63-vector Il{z) = {Qb, z"'). The are taken as projective co- 
ordinates on the complex structure moduli space, with Qi, = dbQ{z). The Kahler potential 
/C for the as well as the axio-dilaton ^ = Cq + ie~'^ is 

/C = - log(z / O A H) - log(-z((^ - ^)) = - log(-zrf ■ E ■ n) - log(-z((^ - ^)) , (2.2) 

J M 

where E is the symplectic matrix E = ^ ^ ^ . The axio-dilaton and complex structure 
moduli take values in the moduli space A^; a correct global description of the moduli space 
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requires that we identify points in Ai related by modular symmetries, as we describe in 
§2.2. 

We now consider nonzero fluxes of the RR and NSNS 3-form field strengths Fs and 
Hs over these three cycles, defining the integer-valued 63-vectors / and h via 

Fs = -{2Tif a' {fa a" + fa+h,,,+i (3a) , H3 = -{27vf a' (ha a" + ha+h,,,+l Pa) . (2.3) 

The fluxes induce a superpotential for the complex structure moduli as well as the axio- 
dilaton ||: 

W= [ GsAniz) = {2n)^a'{f-(Ph)-U{z), (2.4) 

where 6*3 = F3 — (pH^. 

We will be interested exclusively in vacua satisfying the F-flatness conditions: 

D^W = DaW = 0, (2.5) 

where DaW = daW + WdaJC, and we have allowed a to run only over /i2,i inhomogeneous 
coordinates. This is alternately 

(/ - ^h) ■ U{z) = (/ - <l)h) ■ {daU + UdalC) = . (2.6) 

These conditions force the complex structure to align such that the (3, 0) and (1, 2) parts 
of the fluxes vanish, leaving the fluxes "imaginary self-dual," *eGs = iG^- 
The fluxes also induce a contribution to the total D3-brane charge 

In the rest of the paper, we will set (27r)^a' = 1 for convenience. For vacua satisfying 
(p75|), the physical dilaton condition Im (p > implies A^Aux > 0. As the total charge on 
a compact manifold must vanish, sources of negative D3-charge must be present as well. 
For a given IIB orientifold compactiflcation, a flxed amount of negative charge is induced 
by the orientifolds, leading to an effective bound on A^Aux: 

iVflux<i^, (2.8) 

where for instance in a IIB orientifold arising as a limit of a fourfold compactiflcation of 
F-theory [^, L can be computed from the Euler character of the fourfold |^^. Although 
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the number of imaginary-self dual flux vacua is inflnite, the set satisfying (2^) for flxed L 
is in general flnite, as we shall discuss in §2.3. In this paper, we do not restrict ourselves 
to particular values of L associated with a speciflc orientifold, but instead following 
we take L to be an arbitrary positive integer, and focus on developing tools for estimating 
numbers of vacua with certain properties, as a function of the fixed parameter L, in the 
context of this simpler mathematical model. 

In a more realistic model in which one applies the methods developed here to a partic- 
ular compactification of potential phenomeno logical interest, one would need to impose a 
value of L associated with a particular orientifold of the Calabi-Yau. The charge difference 
L — Afflux can be made up by mobile D3-branes; note that antibranes, another potential 
source of negative charge, violate the imaginary self-dual structure. Their inclusion in this 
class of models (which is only consistent with finding a tadpole-free configuration after 
one has stabilized Kahler modes) does not greatly extend the class of fiuxes which must 
be considered, since the number which can be included without inducing a perturbative 



instability is small compared to L p3| , [i4 



Inclusion of Kdhler Moduli 

In general, Kahler moduli are also present in these systems, and they do not acquire a 
potential from the fiux superpotential ( p.4| ). In the leading approximation, however, they 
modify the problem in a very simple way. The supergravity potential takes the form 

y = e^*°*(^^AW^^jW^-3|M^p) , (2.9) 

where i, j run over the complex structure moduli, the dilaton, and the Kahler moduli, and 
ICtot includes the full Kahler potential for all these fields. As described in the Kahler 
contribution to the scalar potential in these models can be shown to be 

g'^'^D^WD^W ^^\W\'^ , (2.10) 

where a, (3 run over Kahler moduli only. The cancellation between the contribution ( |2.10| ) 
and the — 3|VF|^ in ( |2.9|) exemplifies (at leading order in Qs and a') the famous "no-scale 
supergravity" structure [^. As an example, models with a single volume modulus p 



(whose imaginary part controls the Calabi-Yau volume and whose real part is an RR axion 
arising from reduction of the C4 field) have a Kahler potential contribution 

/C(p,p) = -3 1og(-z(p-p)), (2.11) 
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which can be obtained from dimensional reduction along the lines of |5G]; given the p- 
independence of (2^), it can easily be seen that (|2J0D follows. 

In a no-scale supergravity background, the potential becomes positive semidefinite, 
with vacua only arising at y = 0. To obtain a vacuum at ^ = 0, one must solve the 
dilaton and complex structure F-flatness conditions (pT5|). If one has -D<^VF or DaW ^ 0, 
the e^*°* factor in (|2.9|) causes the (positive) potential to fall rapidly to zero as Im p ^ oo, 
leading to rapid decompactification. 

On the other hand, in general one cannot solve D^tW = 0. The equations DaW = 
D(j)W = provide /i^'^(M) + 1 equations in /i^'^(M) + 1 unknowns. Imposing D^tW = 
requires = as well by ( p.lO|) , but since W depends only on cj) and the z°', this 
overdetermines the constraints on these fields, with no solution in general. This means 
that for general choices of the fluxes, supersymmetry is broken by the auxiliary field in the 
p multiplet, while the tree-level vacuum energy vanishes due to the no-scale cancellation 
( p.lOp . Vacua with 1^ = may, however, exist for special fluxes; we discuss these in §2.5 
and count such vacua in various examples in later sections of the paper. 

It is critical to note that the special no-scale structure (|2.10|) is not expected to be 
preserved by quantum corrections. For instance, the first a' correction to the Kahler 
potential |^ ruins the cancellation (|2.10|) . As described in fl^, generic corrections to /C 
and W break the no-scale structure, for example due to a non-perturbatively generated 
superpotential for the Kahler fields, permitting nontrivial solutions to all the equations 
DiW = 0. Therefore, one should expect that under many circumstances, fluxes which 
yield only non-supersymmetric no-scale vacua with D^W ^ 0, yield super symmetric AdS 
vacua with D^W = but W ^ after inclusion of further corrections. These corrections 
are, however, in general difficult to compute exactly. 

Due to this difficulty, we follow the lead of Douglas and collaborators and neglect 
the Kahler moduli entirely, focusing exclusively on the dilaton and complex structure 
moduli from now on. Their philosophy is that the distribution of vacua for the dilaton and 
complex structure moduli alone is more representative of the character of the ultimately 
frozen Kahler moduli then the inevitably broken no-scale structure ( |2.10| ) would be. For the 
questions we are trying to address, concerning the statistics of discrete symmetries which 
act only on the complex structure and dilaton, or the frequency of vacua with W = at 
tree level, the detailed dynamics of the Kahler moduli would be irrelevant in any case. 

This approach should be a reasonable first step in learning about the space of real IIB 
string vacua, as in many cases further corrections will yield AdS or dS vacua starting from 
generic no-scale vacua, along the lines of |l^,|15|,|l^ . 
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2.2. The modular group and vacua 

Obtaining proper statistics of flux vacua requires not overcounting multiple copies 
of the same vacuum that are produced by symmetry transformations. We review these 
symmetries here. 

In the absence of fluxes, a symmetry group Q = S'L(2,Z')^ x F acts on the moduli 
space TVI, where SL{2, Z)ff, is the S-duality of type IIB string theory and F is the modular 
group of the complex structure moduli space. Points on related by Q are considered 
equivalent, and a fundamental domain for the moduli space arises from dividing out by Q. 

For the vacua we consider, the fluxes are affected by Q as well. SL{2, Z)^ acts in the 

ordinary way: given an SL{2, Z) matrix ( have 

+ d' \h J \c d)\h) ■ ^^-^^^ 

Under this transformation (/ — ^ (/ — (/}h)/{c(f) + d), and hence solutions of ( |2^ ) 
are carried into other solutions, while A^Aux ( F^ ) is preserved. The action of SL{2, Z) 
generates a Kahler transformation on W ( p. 4] ) and /C ( p.2|) : 



c d 

a(j) + b f f\ f a b\ f f 



W^AW, /C ^/C- log A -log A, (2.13) 



with in this case A = 1/ {c(p+d) . Kahler transformations (|2.13|) are by deflnition symmetries 



of TV = 1 supergravity; the potential V = e^(|iI>VF|^ — 3|VFp) is manifestly invariant, while 
an appropriate transformation for other flelds, including fermions, is determined by A (see, 
for example, |33]). We shall flnd later that preserved discrete symmetries in various low- 



energy theories also take the form of Kahler transformations. 

Since SL(2, Z) is of inflnite order, we immediately see that the action of this group on 
any given vacuum will produce an inflnite number of copies. To usefully count inequivalent 
vacua, one must mod out by this redundancy. One may do this either by restricting (p to 
a fundamental domain or by placing a restriction on the fluxes; we will flnd both methods 
to be useful. It is worth noting that the matrix —1 is in SL{2, Z), although its action on 
(j) is trivial; hence vacua with the signs of all fluxes flipped are gauge-equivalent. 

The precise nature of the complex structure modular group F will vary from case to 
case. The general form for our examples is that under the transformation of the complex 
structure moduli z"' —>■ z'"", the periods change as 

n(2") A(2") M ■ n(2") , (2.14) 
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where A is a function of the moduh and M is a constant Sp{bs, Z) matrix. ( 2.14| ) then 



induces a Kahler transformation ( |2.13| ), and is consequently a symmetry, provided the 
fluxes transform as 

f^f-M-^, h^h-M-^, (2.15) 

which thanks to the symplectic condition have the same value of A^Aux as if^h), and which 
lead to solutions of ( p.6|) with moduli 4>.,z'"'. Again to count only inequivalent vacua, one 
must mod out by these images under the modular group. We shall have more to say about 
fixing r for the various cases as they arise. 

2.3. Enumerating vacua 

Here we discuss the volumes in flux-space occupied by sets of imaginary self-dual 
vacua given fixed restrictions from the modular group and the D3-brane charge bound L, 
and describe how one can estimate the total growth of vacua with L. We additionally 
review the theory of Ashok, Denef and Douglas ||2^,^ , which takes a statistical approach 
to describing the set of vacua for this class of compactifications. 

The set of integral fluxes fi^hi form a lattice Z^''^ R^^'a Because the F-fiatness 



conditions (|2.5| ) are linear in these fluxes, the set of possible continuous values for the 
fluxes solving these equations and giving rise to dilaton and complex structure moduli 
in any fixed finite size region in the moduli space M. forms a cone C in R^^^; when the 
restriction on the moduli space is to a fundamental domain of Q, we call the cone Cp- The 
cone C intersects the unit sphere S'^^^~^ in a region whose boundary is determined by the 
constraints imposed on the moduli. Because the number of F-flatness conditions is equal 
to the number of moduli, these equations generically have a solution, so that the cone C 
is in general a 263-dimensional subspace. 

We saw in §2.1 that in a given compactification, tadpole constraints require one to 
impose the condition A^Aux < L on the fluxes for some fixed L determined by the global 
properties of the model. Because A^Aux becomes a positive definite quadratic form on the 
fluxes for solutions to the F-flatness conditions (|2.5|), the upper bound on A^Aux has the 
effect of cutting the cone C along a hyperboloid H in R^^^ ; the set of discrete fluxes giving 
rise to moduli in the appropriate region then corresponds to the set of integral points in 
the cone C which also lie inside the bounding hyperboloid H. Generically, the resulting 
set of vacua is finite. 

In particular, if a compact region of finite volume in M. is chosen, then the truncation 
by H of the cone C is a compact region in R^^^ with finite volume. This volume scales as 
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\fL ^ as L is increased, so that for sufficiently large L the number of integral points in 
the volume scales as: 

iVvacua(iVflux < L) ^ L'^^ , (2.16) 

where we have assumed that L is larger than other quantities in the problem, notably that 
L ^ bs- Obtaining formulae like ( |2.16| ) for more specific subclasses of vacua will be one 
of our primary goals; in the future we will simply write A^vacua(-^) to mean a set of vacua 
with Afflux < L. 

Subtleties may arise in this counting, and one of our goals in this paper is to investigate 
examples to help understand under what circumstances the estimate ( p.l6| ) is valid. As an 
example of a case where this asymptotic formula may break down, the boundary of the 
cone Cp in continuous flux space associated with a fundamental domain may intersect the 
codimension one region Z C 5"^^^"^ on which the quadratic form iVflux vanishes. In this 
case the cone truncated by H is no longer compact. As discussed in |2^], in this situation 



the number of lattice points in the truncated cone may not scale according to ( |2.16| ). As 
L is increased, there are lattice points in the truncated cone whose projection onto S'^^^~^ 
approaches Z; the fluxes at these lattice points are generally of order L rather than of order 
\/L. We will encounter an example of this behavior for the "special symmetric torus" in 
§4.2, where restricting to a subset of the allowed fluxes gives rise to an extra factor of log L 
in the number of vacua beyond that which would be indicated from the counting (|2.16| ). 

At this point, we can see why one needs to develop a statistical theory to gain some 
understanding of this class of vacua. Simple orientifold models with L ~ 10*^ and 63 ~ 100 
can be obtained by starting with the Calabi-Yau fourfolds in appendix B.4 of [^^. In 
such models the naive estimate ( f2.16| ) gives rise to ~ 10"^*^*^ vacua. To gain any credible 
understanding of the detailed properties of such a collection, one needs to have a statistical 
description. 

The idea that discrete fluxes would lead to a large number of closely spaced discrete 
vacua was developed in pi|,l2^,p3[ . A precise reflnement of the rough counting summarized 



above, giving rise to such a statistical description, has been proposed in ||2^,^. The 



theory of these authors actually provides much more information than ( |2.16|) , however. In 
addition to giving the scaling of the number of vacua with L, they also provide a formula 
for the distribution of the vacua on the moduli space, under the fundamental assumption 
that L is sufficiently large that the fluxes can be treated as continuous parameters. The 
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most transparent formula is that for an "index" where vacua satisfying (^]^) are counted 
weighted by sgn(det -D^VF): 



-^vacua 



{L)r^L^' I det{-R-u)), (2.17) 
Jm 



where uj and R denote the Kahler and curvature 2- forms on M. , respectively. This formula 
suggests a smooth discretization of flux vacua across the moduli space. For fixed L, how- 
ever, a sufficiently small M.' will begin to probe the "fractal" distribution of vacua coming 
from the true integral nature of the fiuxes, invalidating (|2.17|) . 

Analogous (but more complicated) formulae to ( |2.17| ) are derived in [p5| , |26[] for the 



(unweighted) number of vacua satisfying ( |2.5| ), as well as more general solutions to the 
equations daV = d(j)V = not necessarily satisfying (|2.5|) .i Numerical results consistent 
with the formula (|2.17|) have been obtained in p7| , |29[] , and indications so far are that it is 



correct within its range of validity. 

Our goal will be to compute the numbers and distributions of vacua without making 
the primary assumption that leads directly to the formula (|2.17| ): neglecting the discrete- 
ness of the fiuxes. In part, this will be to allow meaningful independent comparison with 
( p.l7|) . But this is also important because the more detailed questions we are interested 
in — including the relative frequency of vacua with W = 0, and the relative frequency 
of vacua with discrete symmetries — involve overconstrained systems of equations, giving 
new mathematical questions which cannot be answered by the statistical theory. Here, 
we briefiy summarize some useful ways of thinking that emerged in our analysis of these 
models, and may be more generally useful. We then describe the two relevant classes of 
vacua in more detail in subsequent subsections. 

2.4- Algebraic classification of vacua 

One approach to counting fiux vacua is to choose a set of allowed integral fiuxes and 
then to try to solve the equations ( |2.6| ) for each combination of fiuxes. In some very simple 
cases these equations can be solved exactly, as we demonstrate below for the symmetric 



^ These authors count as "supersymmetric" any solutions to ( p.5D , and term these more general 
vacua "non-supersymmetric" . Although the more general vacua would have tadpoles in the no- 
scale approximation after inclusion of Kahler moduli, one can play off these tadpoles against 



further effects to obtain stable, nonsupersymmetric solutions [15|. Due to confusion over whether 



"supersymmetric" means only solving (|2.5D or implies = as well, we avoid the term altogether. 
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torus, but this is in general a difficult task. In simple Calabi-Yau models, for example, the 
periods are given by transcendental functions, and one is often reduced to approximation 
techniques as in ||2^,|9l. While we will also use such techniques here, they cannot easily 
be used to establish the existence of vacua with W precisely vanishing, or vacua at a given 
point in moduli space where a discrete symmetry could be restored. 

Instead of specifying fluxes and then searching for solutions, one can also choose to 
specify a point in A4 and then search for integral fluxes which yield solutions precisely 
at this point. If the points in associated with allowed vacua ("allowed" as defined by 
the interest of the questioner) can be identified to lie in some particular discrete family, 
we can then in principle sum over the allowed vacuum values in A4. We can sometimes 
do this exactly, to arrive at an exact formula for the number of vacua as a discrete sum. 
Sometimes we can proceed by finding generic features of the allowed values of Ai which 
affect the scaling of the number of vacua, and which allow us to estimate the total number 
of vacua in terms of the asymptotic number of points in Ai with those features. We now 
consider this in more detail. 

2.4- 1- Vacuum counting at a point in moduli space: field extensions 

In all cases we examine in this paper, it turns out that the allowed points in moduli 
space where vacua arise are points where (with an appropriate normalization on O) the 
period vector 11(2;) and its first covariant derivatives have components taking values in a 
finite extension of the field Q of rational numbers. For example, one family of extensions 
we encounter on the torus are the imaginary quadratic number fields Q[V~"-D] with D > 0. 
On Calabi-Yau manifolds, we have periods which lie in one of the cyclotomic fields JF(Q)fc, 
which for integer k is the extension of the rationals by powers of the kth root of unity 
a. The algebraic structure of such points in Calabi-Yau moduU spaces has also played an 
important role in the works [|50| , |5T|j5^ , |53| , |5^j55| . 



It seems that the most important characteristic of a vacuum associated with a point 
p E Ai is the degree Vp of the field extension over Q associated with the periods at that 
point. The degree of the extension Q[\^—D] is P = 2, while the degree of the extension 
J^k is given by = (^(/c), where (^(/c) is the Euler totient function,! the number of positive 
integers < k which are relatively prime to k. It can also be computed as 

m = ki[(i-^] (2.18) 



(f>{k) should not be confused with the axio-dilaton 
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where Pi are the prime factors of /c, with each prime factor appearing once in the product. 

The degree of the field extension at a point provides a simple estimate of how many 
flux vacua should arise there, for a fixed bound A^Aux < on the fluxes. The point is that, 
when a single equation in (|2.6| ) fixes a modulus given the fluxes, it imposes precisely one 
condition on the moduli, since they take complex values. If, however, the moduli are fixed 
and one tries to satisfy the equation by constraining the integer- valued fluxes, the number 
of constraints that arise is larger, and is given by the degree V of the extension. 

At a point in moduli space where the field extension JF over Q is of degree P, one can 
view the equations (|2.6| ) (or their generalization including 1^ = 0) as a system of linear 
equations which the integral fluxes must satisfy, with coefficients lying in the field T. Note 
that because these equations are also linear in any vacuum satisfying these equations 
with periods (and derivatives) in JF must also have (p E J-". A general element t of JF can 
be written in an expansion 

t= Y,^^a,, (2.19) 

i=l 

where G Q and cr^ form a basis of elements of JF which are linearly independent over Q. 



For example, for the field extension Q[v — D], -D > 0, we can use ai = 1, (72 = ivD, while 
for the cyclotomic fields JF^, we can take ai = a^^~^, with a a kth root of unity and 
relatively prime to k. We therefore find that solving ( |2.(j| ) at a point in Ai with periods in 
a field extension of degree V can be thought of as solving a system oi rj = T> (62,1 (M) + 1) 
equations in the 263 (M) integral components of the flux vectors /, h. 

This analysis leads directly to a prediction for the generic L scaling of the number of 
solutions associated with any given point in moduli space. Since the quadratic constraint 
-^flux < L implies fi, hi < C(v^), we expect 

iVvacua(i^) ~ L^^^ K = 2bs-V, (2-20) 



as the scaling of solutions to (|2.6| ) at such a point in Ai, and 



iVvacua(i^; W = 0)^ L^/^ 'y = 2bs-r]-V = K-V, (2.21) 

as the scaling of solutions with 1^ = 0. These arguments are only clearly justified for 
generic [W = 0) vacua when k > (7 > 0). This indicates that they will be most useful 
for extensions of very low degree V, though they may apply more generally. 
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In our examples later in the paper, we find some situations in which these formulae 
apply and immediately give a correct estimate of the number of solutions associated with 
a fixed point in moduli space. In many concrete examples, however, subtleties arise. First 
of all, it is not always true that the equations ( |2.6| ) are linearly independent; we will find 
for the example of Landau- Ginzburg vacua of Mg in §5 that the two equations in ( p.6| ) are 
degenerate when the expansion ( p.l9|) of the dilaton is restricted. Secondly, in many cases 
the system is over determined, generically constraining the flux vectors /, h to vanish, such 
as the examples of Landau-Ginzburg vacua M5 and Mio, as well as = vacua for Mg, 
in §5. Thirdly, in some cases one of the constraints on fluxes following from (|2.6|) may 
require, e.g., A^Aux = 0, in which case there can be no nontrivial ISD solution, regardless 
of the counting of degrees of freedom; this occurs for 1^ = vacua in the rigid CY model 
of §3. It is possible other subtleties could arise in other examples. 

2.4.2. Vacuum counting over a region in moduli space: heights 

The analysis we have just described gives an order-of-magnitude estimate for the 
number of flux vacua associated with a particular point in moduli space. In most cases, 
we are interested not just in a single point in moduli space, but in the ensemble of points 
admitting vacuum solutions. Thus, we want to sum the contribution to the total number 
of vacua over a large set of points in A^. In some cases we want to sum over all points 
in moduli space admitting vacua. In other cases, such as on one-parameter Calabi-Yau 
manifolds, we may wish to fix the complex structure moduli and sum over all possible 
values of the axio-dilaton (p. In any of these situations, we require a more refined estimate 
of the number of vacua than (|2.20| ), (|2.21| ). In particular, for two points pi,p2 G 7W we need 
some measure of the relative number of vacua associated with these two points. If these 
points are associated with vacua having periods with the same degree of field extension 
over Q, then even if the estimate ( |2.20[ ) holds, we need to know the relative scale of the 
overall coefficient multiplying L**/^. 

A very useful tool in developing a more precise estimate for the number of vacua at 
some point p ^ Ai is the notion of height, which is a fundamental tool in the modern 
approach to Diophantine geometry (i.e., the study of solutions of systems of equations 
over the integers using algebraic geometry — see for example . For another discussion of 
heights in a somewhat related physical context, see |]51[). The basic concept of a height is 
very simple. It is really just a systematic way of developing an order-of-magnitude estimate 
for problems involving integers and rational numbers. For example, the height of a rational 
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number p/q m. reduced form (so that q have no common divisors) is just the maximum 
of \q\. More generaUy, if we have a n-plet of rational numbers pi/^i, . . . ,Pn/Qn, we can 
think of this as a point 

a:=fl,^,^,...,^')eP"(Q), (2.22) 

V Ql <12 QnJ 



in the n-dimensional projective space over Q. The height of a point ( p.22|) in P"(Q) 



is given by taking a representative of the point where all denominators are cleared and 
common divisors are removed, 

X ~ (mo,mi, . . . ,m,^) , gcd(mo, • • . , m^) = 1 , (2.23) 

and then defining the height of x to be 

H{x) = max(|mo|, • ■ ■ , |?T^n|) ■ (2.24) 

This notion of height is extremely useful in the mathematical study of Diophantine equa- 
tions because it systematizes the notion of order-of-magnitude in a way which is natural 
to number-theoretic problems. While this may seem rather arcane, we will only use this 
notion in a fairly straightforward physicist's sense, in order to compute the scaling behavior 
of simple systems of equations. 

The basic idea is the following. Imagine we have a linear equation 

Ao + + ^A2 + ... ^Ar, = , (2.25) 

for fixed values of the rational coefficients Pi/qi, and we wish to estimate the number of 
solutions of this equation for integers Ai with \Ai\ < N, for large (that is, the height 
of the solution (Aq, . . . A^) is bounded at A^). As above, the scaling is clearly generically 
as A^", but we are also interested in the overall coefficient. We can understand this by 
considering the height of the coefficients, treating them as a point ( |2.22|) . 



Clearing denominators as above, the solutions to the equation (|2.25| ) with rational 
coefficients are clearly the same as the solutions of the equation with integral coefficients 

moAo + miAi -\ ninAn . (2.26) 

It is easy to see that generically the number of solutions of this linear equation goes as 
N'^ /H{x) where the height H{x) of the coefficients is defined through ( p.24|) . (For example. 
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assume that mo is the largest of the integers m^. Then the remaining sum over rriiAi, i ^ 
is of order HN or less, and generically 1/H oi the possible values of that sum will be 
divisible by mo; thus, we expect of order N^/H solutions.) 

We can generalize the notion of height slightly to include the field extensions over Q 
discussed above. If we have, for example, a modulus ^ G JF, where JF is a field extension 
of degree T> over Q, we can as in ( p.l9| ) consider ^ to be a P-plet of rational numbers. 
Considering these numbers to define coordinates in rational projective space, we can define 
the height H{(p) as above. Similarly, we can define the height of a polynomial 

Gnx'^ ^ aix + ao G Q[x] , (2.27) 

with rational coefficients = Pi/qi to be the height of the corresponding point 
(ao, . . . , ttn) G P'^(Q). We will use this generalization of the notion of height in particular 
in the example discussed in Section §4.3. 

One helpful property of heights is that they behave well under multiplication. If 
a, (3 & J-" have heights H{a), H{(3) then generically 

H{a(3) ~ H{a)H{l3) . (2.28) 

Similarly, if polynomials f{x),g{x) G Q[x] have heights H{f),H{g) then generically 
H{fg) ~ H{f)H{g). This statement can be made more precise (see for example []5B|), 
but we will simply treat it as an order-of-magnitude estimate. 

We will also find it very useful to estimate the number N{n, H) of points x G P"^(Q) 
of height < H as H^^^ . Since generically ti + 1 integers of order H have no common 
divisor, this estimate follows directly from ( p.23| ). A more precise calculation shows 
that 

on 

iV(n, H) = —— + 0{H-) , (2.29) 
C(n + 1) 

where the sub leading term goes as H log H for n = 1. The factor of 1/C(n- + 1) = 
Ylp{l — l/p""*"^) where p runs over all prime numbers serves to remove the fraction of 
integers with common divisors — for example, 1/2"^+^ of all (n + l)-plets of integers share 
the divisor 2, 1/3"^"'"^ of all (n + l)-plets of integers share the divisor 3, etc. 

We can now describe a general estimate for the number of vacua where we sum over 
points in the moduli space for the axio-dilaton and/or the complex structure moduli. 
Assume, for example, that we fix the complex structure moduli z so that the periods take 
values in a particular field extension JF^ of degree V over Q, and we consider all values of 
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the axio-dilaton lying in Tz- For simplicity, let us take the periods to be of order 1. For 
each possible value of the axio-dilaton, we have as above r\ = 1^(62,1 + 1) linear equations 
which must be satisfied. The total number of fluxes is 263, with each flux of order \fL. If 
the axio-dilaton has height if, then each linear equation is solved by a fraction of order 
1/{H\/Li) of the possible fluxes. Generically then, the number of vacua at some fixed (j) of 
height H[(f)) is given by 

N^^,UL;(t),z)^ . (2.30) 

One subtlety that will be important to us is the following. For each equation of the form 
f ■ li = (j)h ■ where 11 is the period vector or one of its covariant derivatives projected 
onto a basis element of JF^, generically h ■ 11 will be nonzero, and the height of (j) will 
constrain the possible solutions for / by leading to the counting in (|2.3C1| ). However, 

it is also possible for both sides of the equation to vanish separately; this is supressed by 
a factor but enhanced by if, since the height of (j) provides no constraint. In most 

cases the generic solutions will dominate, but for H > \[Tj solutions with / ■ 11 = /i ■ 11 = 
may contribute to, or even dominate the generic solutions enumerated by ( p.30|) . As we 
discuss in the next subsection, 1^ = vacua involve precisely these special solutions for 
the DaW = equation. 



We can also sum (|2.30| ) over all values of the dilaton. Since there are on the order of 
allowed values for the axio-dilaton of height if, the total number of vacua for all (p 
would be 

iVvacua(i^; Z)^ Yl H^'"' L^''''/^ , (2.31) 
H=l 

where the upper bound ifmax is the maximum dilaton height consistent with our constraint 
on the fluxes. (In all cases we are studying in this paper, either the degree of the extension 
is 2, in which case ifmax = L, or the sum converges so rapidly that the precise upper bound 
is irrelevant.) 

Summing over various complex structure moduli z is more complicated, particularly 
when we wish to sum over a range of z involving different field extensions. Understanding 
the proper techniques for performing this type of sum is an important question for the 
future. 

These kinds of arguments can be generalized, for example to solving for vacua with 
= 0, which we do in the next subsection. If the preceding discussion seems a bit 
abstract, the reader may find it useful to go through the example discussed in section 
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§3.3, which gives a concrete example of how to use this kind of analysis. We will use the 
considerations described in this subsection as one tool for estimating the number of vacua 
in many of the simple models we study. 

In many cases subtleties arise that must be taken into account. In particular, the 
formulae ( |2.30|) , ( |2.31|) were derived under the assumption that the coefficients in the 
linear equations satisfied by the fluxes are not correlated in any way. In many cases this 
is not a valid assumption — for example the periods may satisfy algebraic relations which 
need to be incorporated in a correct analysis of the space of solutions. Furthermore, in 
some cases logarithmic factors arise which correct ( p.31|) . Thus, these formulae are not 



applicable in all models. Even when these precise formulae are not correct, however, the 
considerations we have used in deriving them provide a useful paradigm for framing the 
enumeration problems of interest. By understanding precisely how these formulae are 
violated in specific models, we can often develop a more concrete analysis giving precise 
answers. 

2. 5. Vacua with = 

For several reasons, it is of interest to count vacua with = in the no-scale ap- 
proximation. While we will reduce this to a clean mathematics problem in the models and 
approximations we use, it may be useful to provide the underlying physical motivation(s) 
to study these vacua. We will give two motivations, and then describe why the problem 
cannot be answered by the statistical theory of §2.3. 

1) By standard nonrenormalization theorems, one expects that the condition = 
persists to all orders in perturbation theory, with possible corrections only being generated 
by non-perturbative effects in the a' or Qs expansion. Such effects are exponentially small in 
parameters (like bare gauge couplings) that one may expect to have a uniform distribution. 
Hence, the distribution of W values, in any vacua which have Wtree = 0, will quite plausibly 
satisfy 

iVvacuadW^lV^p < e) ~ , (2-32) 

log(e) 



for given e. This in particular makes small values of W natural. Given the formula 
for the vacuum energy in A/" = 1 supergravity, one sees that any SUSY-breaking F-terms 
in a vacuum with small cosmological constant A need to satisfy 

A«A4^L^~L_^. (2.33) 
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In a scenario where the entire W ^ is generated by nonperturbative dynamics, it is very 
plausible that the same holds true for any SUSY breaking F expectation values. In this 
case, one could also expect the Fa and VEVs to be distributed according to a formula 
analogous to ( p.32|) . 



It was argued in |]3^ that vacua with the distribution ( |2.32| ) for W and the SUSY 
breaking F-terms could give a significant or even dominant contribution to the total number 
of vacua with A ^ Mp, as long as the conditions leading to (|2.32|) (namely, W = ai 



tree-level, together with suitable nonperturbative effects) are not too rare. Since this is 



contrary to the intuition coming from other arguments [p2| , |32|j2^ , |3l| which point to high- 
scale SUSY breaking as being preferred in a landscape picture, finding the fraction of 
W = vacua is rather well motivated. 



2) It follows from the work of that there is a deep connection between the existence 
of R-symmetries and the possibility of spontaneously breaking supersymmetry in the true 
ground state of a supersymmetric field theory. Assuming the superpotential is a generic 
function of the fields consistent with all symmetries, the existence of an R-symmetry is 
a necessary condition for SUSY breaking, and a spontaneously broken R-symmetry is a 
sufficient condition. Since unbroken R-symmetries imply W = 0, the study of W = 
vacua which have W generated by nonperturbative effects in accord with ( |2.32| ) , is closely 



related to the problem of studying R-symmetries in string models ||33|| . We note, however, 
that it is well known that the superpotential is often not a generic function of the fields 
consistent with all symmetries in supersymmetric field theories, so the correspondence is 
not precise. 

Having motivated their study, consider now the equations describing aW = vacuum. 
As we have described, imposing W = as well as the F-flatness conditions D^pW = DaW = 
imposes /i^'^ -|- 2 conditions on the /i^'^ -|- 1 moduli; in general, this system cannot be 
solved. We also have, however, 263 (M) discrete variables at our disposal - the choices of 
fi and hi. At least at special loci in moduli space in special models, the freedom obtained 
by varying the integer flux choices does suffice to give solutions to the system p,^PJ^,|[. 

Hence the existence of W = vacua will depend crucially on the choice of fluxes. For 
continuous fluxes, the system would in fact not be over determined - one could view the 
fluxes as additional fields, and roughly speaking solve the extra equation by solving for one 
of the fluxes in terms of the others. This is not possible, however, in the real string theory, 
where fluxes are quantized. The mathematical problem of determining which fluxes are 
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associated with a vacuum with this further constraint becomes one of solving a system 
of equations over the integers. GeometricaUy, we are considering the intersection of the 
cone C discussed in §2.3 with the hypersurface = 0. How many, if any, integral points 
lie on this hypersurface is a problem which cannot be solved using a simple continuous 
approximation to the fluxes. But the continuous approximation is a necessary step in 



deriving the statistical theory of [|2^,^. Therefore, we do not expect this theory to be 
easily generalized to capture properties of the special subspace of = vacua. 

Instead, we must use more number-theoretic methods, as outlined in the previous 
subsection. The most basic question is, what fraction of the solutions have = 0? 
Generalizing the results of §2.4, one finds that a subtlety arises: it is easy to see that the 
equations W = D^j^W = together are equivalent to the axio-dilaton-independent relations 
/ ■ n = /i ■ n = 0. Consequently instead of (|2.30|) one finds 

^63_(^+X))/2 

iVvacua(i^; (P.z;W = 0)r^ , (2.34) 

where there are more powers of 1/VL due to the W = equation, but actually fewer 
powers of l/H, and thus 

-^max 

iVvacua(i^; z;W = 0)r^Yl H^^-^ 1'^-^^+^^' . (2.35) 

We will look for W = vacua and determine how many of these vacua there are in several 
examples. In the concluding section, we will speculate about the answer in more general 
models, informed by our present results. 

2.6. Vacua with discrete symmetry 

The other special kind of flux vacuum we are interested in is one with an enhanced 
discrete symmetry. In this subsection we discuss these vacua, beginning by describing 
precisely what is meant by the statement that a vacuum possesses an enhanced symmetry. 

Each choice of fluxes deflnes not just a vacuum value for the moduli flelds, but a 
complete low-energy theory for these flelds, with the fluxes appearing as coupling constants. 
A proper symmetry of this low-energy theory must be a transformation of the flelds, but 
not the couplings, which are not treated as dynamical quantities. 

Hence, when we speak of "a vacuum with an enhanced symmetry" this is shorthand for 
a set of fluxes for which the low-energy theory of the moduli, including the superpotential. 
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is symmetric under a transformation of the moduli fields alone, with the fluxes remaining 
inert. Once a symmetry is established in the low-energy Lagrangian, one may ask whether 
the vacuum values of the moduli fields spontaneously break it. 

We now discuss the framework within which all the symmetries we discuss emerge. The 
full string theory has a number of symmetries acting both on the fields and on the fiuxes 
as well. Let us consider a group G of such symmetries; a subgroup of the modular group 
G d Q will be our primary example. In general G-transformations cannot be considered 
low-energy symmetries since they act on the fiuxes. One can, however, associate with G 
the set of transformations H that act on moduli just as G does, but do not change the 
fiuxes. In general H will not be a symmetry of the system. It is natural to ask, however, 
whether there are circumstances under which H is indeed a good symmetry. 

Since the kinetic terms for z"" are independent of the fiuxes, G-invariance implies 
they are invariant under H as well. Hence the question of the validity of H reduces to a 
question about the potential Vf^h{4>^ associated with the fiuxes / and h. Let /, /i, ^, 
be mapped to /', /i', (j)' , z""' under G and to /, h, (f)' , z""' under H. i7-invariance requires 

yj,;,((^',^«') = ^M(<^,^")- (2.36) 

The assumption of G-invariance, however, means Vf'^h'{(p'7 ^"') = z"'): which implies 

the condition ( |2.36| ) is equivalent to 

Vf>,h'{(l>.z^) = VfM(l>,z''). (2.37) 

That is, the condition (|2.36| ) that the potential is invariant under the if-transformation, is 
equivalent to a condition that the potential is equal as a function of the moduli to another 
potential generated by G-image fiuxes. 

One condition for two potentials to be equal, of course, is that they have the same 
minimum. As discussed in the previous section, the DW = vacuum is the unique 
absolute minimum of the no-scale potential; hence fiuxes /, h and their images /', h' can 
only satisfy ( p.37|) if they lead to the same solutions to ( ^75|) . But by acting with G on 



we know that the vacuum solution (po, Zq with fiuxes /, h is mapped to a vacuum solution 
(Po'^Zq' for the fiuxes f',h'. Thus we see that the potential can only be iif-invariant if 
't'o = 't'o': = Zq\ that is, if its vacuum sits at a fixed point of G on the moduli space. 

All our examples of discrete symmetries will be of this type: if the fiuxes are such that 
the vacuum value of the moduli lies at a fixed point of G, than the low-energy theory may 
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be invariant under the corresponding transformation H of the moduh alone. For example, 
for G d Q there are several candidate fixed points. Fixed points of SL(2, Z) preserving 
G = Z2 and G = Z3 occur on the dilaton moduli space; the Landau- Ginzburg point in 
the one-parameter Calabi-Yau models Mk we consider in §5 is a fixed point stabilizing a 
G = Zfc C r. We shall find a number of vacua sitting at these fixed points that have an 
associated H symmetry, though we shall also learn that this does not always happen; the 
vacuum can sit at the fixed point of G without the H symmetry existing. 

These vacua can be counted using the techniques of §2.4. In particular, since they 
occur for fixed values of the moduli, this analysis is well-suited. For the cases where they 
sit at a particular point in the dilaton moduli space, instead of summing over the heights 
to obtain (|2.31| ), one instead evaluates (|2.30|) at the relevant value of H{(f)o). We shall have 



several examples of counting vacua with enhanced symmetries. 

One should note that the argument given above implicitly assumes that the vacuum 
solution Zq is unique. If there is a fiat direction at DW = for fixed fiuxes, it is possible 
to satisfy ( |2.37|) with {(po, Zq) 7^ {(po', Zq'), as long as both points lie along the flat direction. 
In the usual case that the vacua sit at fixed points of G, the preserved symmetries arising 
in these examples are never spontaneously broken, but for the case of a fiat direction they 
may be. 

Besides the modular group Q, we shall mention one other source of potential symme- 
tries, namely the complex conjugation transformation C. Given a solution /, h, (po, Zq to 
( P^) , one may take the complex conjugate to generate a new solution, provided one can 
write 

U{T) = U ■ n(±^") , (2.38) 
for some matrix U and choice of sign for ±z". One then has a second solution to ( p.6| ): 

f = f . U ^h' = -h ■ U ,(t>' = , z"" = ±T , (2.39) 

where we arranged the signs on h' and (p' to ensure that Im cp' > 0. The transformation 
( p.39|) is an antiholomorphic involution of the space, and sends W{<p, z*) — * W{(j),z^). Such 
a complex conjugation transformation will appear for all of our examples. Note that unlike 
the modular group C relates inequivalent vacua, and one should not mod out by it. 

Symmetries of the low-energy theory can also descend from C, as we demonstrate in 
§3.5.2, §4.3.4 and §5.2; we will find that in general these transformations of the moduli 
alone act as combinations ofW^W and a Kahler transformation. As is familar from the 
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analogous situation in the heterotic theory [^|, such a transformation is only a symmetry 
of the low-energy action when combined with spacetime parity. In the heterotic case, chiral 
matter was conjugated under an antiholomorphic transformation, meaning the symmetry 
corresponded to spacetime CP; we expect a similar identification here, although we do not 
consider matter explicitly. 

We shall find in some cases a relation between the = vacua and discrete sym- 
metries: we shall in some cases (but not always) be able to "explain" = such vacua 
in terms of a preserved discrete R-symmetry under which VF(^o) We will 

comment on these symmetries when they arise. 

We have argued that symmetries of the complete theory acting on fiuxes and moduli 
may restrict to symmetries of the moduli alone when the fiuxes produce a fixed point 
vacuum. It is possible to imagine other symmetries of the low-energy theories of the 
moduli that are not restrictions of this type, but we have no such examples. Certainly 
transformations associated to the modular group and to complex conjugation are by far 
the most natural, and we will accordingly focus on them. 



3. Flux vacua in a rigid Calabi-Yau 

The simplest model in which to study fiux vacua is a rigid Calabi-Yau with no com- 
plex structure moduli. This model was analyzed in and the scaling ( ^.16| ) was 



reproduced. We review this model in detail here, adding comments on the distribution of 
vacua as well as analysis of vacua with = and with discrete symmetries. Most ques- 
tions about this model can be answered fairly easily analytically, so that it is a good trial 
laboratory for methods which can be generalized to study more complicated models. We 
analyze the number and distribution of vacua in this model using the methods described 
in §2.4, giving a simple example of the application of these techniques. For comparison, 
we also review the analysis of [25, ^| using the continuous approximation method. While 



this model admits no VF = vacua, it does admit vacua with discrete symmetries. We find 
several classes of such symmetric vacua, with Z2 and complex conjugation symmetries. 
We enumerate the vacua in each class. In the case of complex conjugation, the enumera- 
tion of vacua with this symmetry requires a number-theoretic analysis going beyond the 



continuous methods of 25 
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3.1. Superpotential and equations of motion 

For a rigid Calabi-Yau, there are no complex structure moduli, so 63 = 2 and a 
symplectic basis for H^{M) consists of a single A cycle and a single B cycle. We take the 
periods of the holomorphic three-form Vt to be 

[ n = i, [ n = i. (3.1) 

Jb J a 

The resulting flux superpotential is 

W = A(I> + B, A = -hi-ih2, B = h + ih. (3.2) 

The D3-brane charge associated with these fluxes is 

iVflux = hh2- hif2 , (3.3) 
and the equation for a vacuum in this model has the form 

D^W = A'$ + B = 0. (3.4) 
Solving for given the fluxes, one finds 

0=-|. (3.5) 
The requirement that Im > directly implies ATgux > 0. 

3.2. Generic vacua 

For the rigid model the entire modular group is ^ = SL{2, Z)^. To avoid overcounting 
SL{2, Z) copies of vacua, we wish to gauge-fix Q. We can fix the gauge either by choosing 
a fundamental region for cp or by making a canonical choice for the fluxes. We will use 
both of these approaches to gauge fixing for the various models studied throughout the 
paper, so we review both here. 

The first approach to gauge fixing is to choose (p in the usual fundamental region of 
SL{2,Z), 

(j)eJ^D = {z: -1/2 < Re ^ < 1/2, \z\ > 1, \z\ 7^ 1 for Re ^ < 0}. (3.6) 

This does not entirely fix SL{2,Z); one must remember that fluxes with opposite sign 
(which sit at the same point in .Fd) are SL{2, Z)-equivalent and hence there is an overall 



23 



weighting of 1/2 to properly count vacua. In addition, there are two points on JF^ that 
stabihze a subgroup of SL{2,Z), and hence one must further mod out by a factor 1/2 
for vacua at the Z2 fixed point (p = i, and by 1/3 for vacua at the Z3 fixed point (p = 
exp(z7r/3). (These fixed points have the potential to support discrete symmetries, as we 
will discuss in §3.5.) As an example of correct counting of flux vacua using these weights. 



at Afflux = 1 there are four flux combinations leading to solutions of (|3.4| ), {A,B) = 
(1, i), (z, —1), (—1, — z), (— i, 1). Together these represent a single physical vacuum aX (p = i. 

The other way to fix SL{2, Z) is to constrain the fluxes to take a canonical form. For 
example, we can require (as in [^) that 

/ii = 0, < /2 < /i2 . (3.7) 

The tadpole condition Aflux < L then gives /i/i2 < L. 

It is easier to count the total number of flux vacua using this second representation, 
in which each flux combination satisfying the constraints represents a single vacuum. The 



total number of inequivalent vacua with Agux < in the fluxes is then [^5| 



L L ^ T 2 

A_„a(L) = E = E ^(^) = E^L^J ~ i^^' • (3.8) 

m=l k\m m=l k=l 

where a{m) is the sum of the divisors of m and \_x\ denotes the greatest integer < x. 
3. 3. Distribution of vacua 



It was shown in |2^ that when L is sufficiently large that the fluxes can be considered 
as sampled uniformly from the continuum, the density of flux vacua on the moduli space 
of (p goes as ( |2.17| ) L^/(Im cp)"^. It was furthermore noted in that the distribution of 



vacua at flxed L features "voids" containing no vacua except for an accumulation at their 
centers; for example, the points (p = ki, k & Z are the centers of the largest voids. We now 
describe these results using three different approaches. While aside from some small new 
observations regarding the void structure the results we describe here are not new, and 
essentially recapitulate what was found in [25,26], this simple model provides an excellent 



framework in which to compare the different methods we use throughout the paper in more 
complicated examples. 
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Scaling analysis 

First, let us consider the scaling analysis discussed in §2.4. We can directly apply the 
arguments leading up to ( |2.31| ). In this case our field extension is JF = Q[z], the complex 
rational numbers, with degree V = 2. We have /i2,i = 0, 63 = 2, so the 4 fiuxes satisfy 
rj = 2 equations. Given a value (j) = [n -\- im) / k for the axio-dilaton, which has height 
H{(p) = max(m, n, k), we have from (|2.30| ) 

iV,,,,,(L,(/>)~-^^. (3.9) 

This gives an estimate for the number of vacua at a fixed axio-dilaton (p, as a function of 
the height of (p. To get the total number of vacua, we must sum over (p. There are order 
values of (p with height < H, and therefore order with height of order H. This 
gives the estimate from ( p.31| ) 

L r 

A^vacua(L)~5^if2— -^~L2. (3.10) 

This estimate gives us information about the scaling of the total number of vacua as a 
function of L, and furthermore tells us something about the distribution of values of (p. 
This distribution is expressed here in terms of the height function H, however, rather than 
the more familiar quantity Im cp. Unlike Im cp, which is a smooth continuous function on 
the axio-dilaton moduli space, the height function is a number-theoretic object which only 
is defined at special values of (p and which cannot be continued to a continuous function. 
To compute the distribution on moduli space as a function of Im we need a slightly 
more detailed analysis. This can be done fairly easily as follows. 

For Im (p = m/k ^ 1, in the fundamental domain the height of (n + im)/k is just m. 
The number of values of (p of height H with imaginary part m/k > y then goes as (H/y)'^, 
so the density of (p at height H goes as H^/(lm. (p)^. Each term in our linear equations is 
now of order H^jL/lm. (p, since B /A is of order Im (p, so the expected number of vacua at 
fixed Im cp goes as 



L 



iVvacua ~ 2^ ^^3 ^^^^2 . (3.11) 

Thus, we see that the methods described in §2.4 allow us to compute both the total number 
of vacua and their distribution, up to overall constants. To determine these constants, a 
more careful analysis must be done in which an explicit gauge-fixing is carried out. We 



25 



now compare two more detailed methods of analysis using the two different gauge fixings 
discussed above. 



Continuous analysis 



For comparison, we review the method used in to analyze this model by approxi- 
mating the fluxes as continuous variables. Consider the gauge fixing (p G JF^. The number 
of complex numbers z = Zi + iz2 with -21,-22 ^ Z of magnitude \z\ goes as 27r|-2;|(i|-2|. For 
large \z\, the phases of these numbers are uniformly distributed on the circle. Thus, we can 
estimate the number of pairs A, B (|3.2|) giving rise to any particular cj) in the fundamental 



domain as follows. Taking A = |A|e*^, B = \B\e'^^^^'^^ — e^^ j3 we have 



L > \A\\B\ sin?/;, Im 



sini/', (3.12) 



so \A\ < y^L/lm (j). Thus, the distribution of vacua scales as 

1 r\^L/lm 4) I- 17-2 

iVvacua(L;</>) ~ -y^ 27r|A| d\A\ j (3 5\(t^ - —(3) = 4^^^^, (3-13) 

where the overall factor of 1/2 arises due to the S = —1 symmetry mentioned above. The 
integral over moduli space gives 7r^L^/12, in agreement with the computation ( p.8|) . 

We can also understand the voids easily from this point of view. For any A and y G Z 
we can choose B = iyA as long as |Ap < L/y. Thus, we expect on the order of nL/2y 
vacua at each integral point (p = yi. These solutions appear for any A. The nearest other 
vacua must have a B which differs by at least one, and which is thus separated by at least 
a distance of \fyjL. This simple analysis thus gives us an easy way of deriving the density 
of vacua on moduli space as well as the void structure. 

Discrete analysis 

We now review a more discrete approach to the analysis, using the second approach to 
gauge fixing, where we fix the fiuxes as in (|3.7[) . This approach was discussed in |j2^, but 
is closer in spirit to the height-based analysis above than the continuous approximation 
approach. With the gauge fixing (^77|), the moduli (j) associated with the given fiuxes fit a 
very simple pattern, lying on the lattice of points 



/2 + ^/l 
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(3.14) 



with < /2 < fih2 < L. By using the SL(2,Z) symmetry — > — 1, we can shift 
the set of aUowed fluxes to he in the band —1/2 < Re (p < 1/2. While the total number 
of resulting flux solutions in this gauge fixing is easy to compute and is given by (|3.8|), 
the distribution of fluxes is less obvious. The fluxes just described can be thought of as a 
set of "layers" of h2[L/h2\ vacua, each with density /i2, arranged uniformly on the region 
l//i2 < Im (j) < [L/h2\/h2. This would seem to give rise to a distribution in Im (f) scaling 
as 

h2 = l ^ 

Note, however, that in this gauge fixing most of the vacua lie in the region < 1. These 
points must be moved back into the fundamental domain with an 5'L(2, Z) transformation 
in order to compute the real density. As the vacua in the flux gauge fixing are most dense 
in the region with small Im (p, these vacua lie in many different modular regions. If the 
distributions of vacua on the different modular regions are essentially uncorrelated, the 
image on the fundamental region should give rise to a distribution which has a functional 
form invariant under modular transformations — this is presumably one way of explaining 
the l/(Im (j))'^ dependence of the distribution in this gauge fixing, but this explanation is 
certainly less transparent than in the direct gauge fixing of (p in the fundamental region. 

We can also consider the distribution of voids in the fiux gauge fixing. For the layers 
described above, the appearance of the voids is quite transparent: each layer contains a 
point at each integral point (p = yi, and adjacent points on that layer are separated by a 



distance l//i2 with < L, so that the size of the hole is clearly again y/Im (p/L. The 
situation is complicated by the images of the vacua in the region \(p\ < 1. Note, however 
that the Jacobian of the map from this region into the fundamental domain is greater than 
one, and therefore increases the spacing of the vacua. 

Comments on exact analysis vs. asymptotics 

In the preceding discussion we have given analytic estimates for the rate of growth of 
the total number of vacua and of the distribution of the vacua on moduli space. These 
analytic estimates are based on two types of analysis. In the continuous analysis approach 
of |2^, the estimates are good only when we have a system of equations which is not 
over constrained, and when the fiuxes are large enough that their integral discreteness can 
be ignored and we can use a continuous approximation. In the discrete type of analysis, 
the continuous limit is based on asymptotic formulae for number-theoretic properties of 
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the integers, which can be encoded in functions hke a{k). In this model, we found that 
both approximation methods lead to the same asymptotic formula for the number of vacua, 
( p.l3|) . An important question is how large the fluxes must be for these approximations 
to be valid. Some discussion of this issue was given in |]25|J2^. For the rigid CY model 



discussed here, it is easy to put the exact discrete equations on a computer and demonstrate 
that the approximations are already valid at quite small values of the fluxes. 

It is important to note, however, that in general the discreteness of the fluxes will 
strongly affect any equations involving quantities which are not summed over all integers 
below some bound L, but rather involve speciflc properties of a particular integer A^. For 
example, as discussed above the total number of flux vacua in the rigid CY model at a 
particular value of Afiux is given by 

iVvacua(A^flux) = ^(A^flux) = k. (3.16) 

I Afflux 

This function of A^Aux is not smooth and depends on the number theoretic properties of the 
argument. Both this function and its fluctuations are of order Afiux ^ this is clearly illus- 
trated in the flgure below. On the other hand, the total number of vacua with A^Aux ^ L 
is much better behaved-this function, as described above, scales as L^, with fluctuations 
which decrease in relative scale as L ^ oo. Thus, in comparing precise analytic or numer- 
ical calculations to asymptotic estimates, we will be much better off when we compare to 
quantities like the total number of vacua with A^Aux < L, involving inequalities, rather than 
quantities like Avacua( Afflux), which are much more dependent on factorization properties 
of their arguments. 
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Fig. 1: Plot of A^vacua (Afflux) Up tO Nfiu. = 1000. 
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3.4- Vacua with W = 

In this simple model it is easy to show that the system does not admit any W = 



vacua. W = implies, using (|3.2| ), that Acp + 5 = 0, which according to the analysis of 
§2.4 should constitute X'(Q[i]) = 2 additional constraints on the four fluxes, which naively 
should lead to solutions. It is easy to see, however, that combined with DW = 0, these 
particular constraints force A^Aux = 0, leaving no physically acceptable flux vacua. This 
is an example of a subtlety in the counting of constraints, where naively the number of 
solutions should be nontrivial but the special form of the constraints removes them all. 
One can further note that satisfying W = DW = requires Im cp = 0, which is not an 
admissible point in the moduli space. 

3.5. Discrete symmetries 

Following the arguments of §2.6, we look for enhanced symmetries in theories whose 
vacua sit at fixed points in the axio-dilaton moduli space. For subgroups of SL{2,Z), 
there are fixed points preserving Z2 and Z3; in the fundamental domain these are located 
dX (f) = i and (/> = exp(7rz/3) respectively. Fixed points of complex conjugation C, which in 
this case just fiips the sign of the RR scalar, lie on the imaginary (j) axis. We consider each 
of these types of discrete symmetries in turn. 

3.5.1 Z2 C SL{2,Z) 

For the rigid CY model, fiux vacua exist precisely at the Z2 fixed point (p = i. This 
point is invariant under the transformation S = ( ^ J . As S"^ = — 1, the so-called Z2 is 



^-1 0^ 

really a Z4 when action on the fiuxes is included. As discussed above, therefore, the vacua 
at (p = i will come in quartets, each set of four standing for one inequivalent vacuum. 

The point (/> = z is an element of Q[i] of height 1. From the scaling analysis of §2.4 we 
expect to have on the order of L vacua at this point. More specifically, fiuxes compatible 
with (p = i are of the form 

^i = -/2, /i2 = /i, iVfl,, = /f + /|. (3.17) 

Dividing out by Z4, we find that the number of vacua with A^Aux < L at (p = i is then 

iVvacua(i^,^ = ^) -TT^ . (3.18) 

The ratio of vacua at the Z2 point to generic vacua is thus asymptotically 3/{nL). 
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We next ask whether the low-energy theory inherits the fixed point symmetry. As 
mentioned in the previous section, this would mean the potential is invariant under the 
restriction of the fixed point symmetry transformation to (f) alone, with fiuxes invariant. 
The SL{2, Z) transformation —1 is trivial acting on (p, and hence we are left asking about 
H = Z2 on (f) generated by S: 

(j)^-]-. (3.19) 



It is easy to see that ( |3.17p is equivalent to B = iA. Consequently for these fiuxes the 



superpotential is simply 

W{(j)) = {(l) + i)A. (3.20) 
It is then straightforward that the transformation ( |3.19| ) gives 

W{(p) ^ ^W{(p) . (3.21) 



Hence up to a phase, W{(p) transforms in the canonical way (|2.13|) under SL{2, Z), implying 
that the low-energy theory is indeed invariant under Z2 for any fiuxes leading to a vacuum 
at (f) = i. 

This is our first example of a preserved discrete symmetry: all vacua sitting at the 
fixed point "inherit" a Z2 symmetry acting on the moduli alone. As we shall see, for more 
complicated examples, this will not always be the case. 



There is an alternate derivation of ( 3.21| ) that will generalize to more complicated 



models. We note that another way of expressing the constraint (|3.17D is h = f ■ S, where 



the matrix S is acting now acting on the space of periods, i.e. it rotates fi into /2 rather 
than fi into hi. One can then show that under ( p.l9|) , 

{f-<ph)-.^{f-<Ph)-S, (3.22) 

Using SIl = in, (|3.21J ) then follows. This kind of argument will be useful later when the 
periods are more complicated, but are known to be eigenvectors of certain transformations, 
namely the modular transformations F. 

Note in ( |3.21| ) that because in the vacuum (f)o = i, we have W{(f)o) ^ W{(j)o)- Hence 
even though this Z2 acts on the superpotential it does not constrain W to vanish in the 
vacuum; as noted in the previous subsection, there are no such W = vacua in this model. 

As discussed previously, the Z2 point is surrounded by a "void", whose size shrinks 
as as L is increased. So at reasonably small values of L, there are no vacua nearby 

with an only slightly broken Z2. 
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There is also a fixed point on dilaton moduU space preserving Z3, but it is easy to 



show that there are no vacua precisely at the Z3 point, simply because (j) in (|3.5|) has 
rational real and imaginary parts. In the language of §2.4, this value of (f) does not lie in 
the field extension JF over the rationals Q defined by the periods. We shall find later a 
model with one complex structure modulus where the Z3 point on the (/)-plane becomes 
accessible. 

If, however, one is interested in only approximate symmetries, one may consider vacua 
within some small neighborhood of size e around the Z3 point; the density of vacua there 
will scale like e^L^. For reasonable values L and a small symmetry breaking scale e, the 
Z2 point captures many more vacua than the Z3 point. Note that as soon as e > 
one starts including vacua from outside the void around the Z2 point. So at any fixed e, 
for very large L there will be approximately as many Z2 vacua as Z3 vacua. For smaller 
values of L so e < l/\fL (but still L ^ 1, where the continuous approximation to the 
fluxes should apply), the Z2 vacua are more numerous. 

3.5.2 Complex conjugation Z2 

Finally, we consider the Z2 associated to complex conjugation C. For the rigid CY 
model ( |2.38| ) is satisfied somewhat trivially with C7 = CT3, the third Pauli matrix, so one 
can map one vacuum to another (inequivalent) one via 

fi = fi , f2 = -f2 , hi' = -hi , h2' = /i2 , (p' = -'(p. (3.23) 

The fixed points on the (p-plane are of course aX (j) = i(3, which are vacua for fiuxes 

fi=Ph2,f2 = -Phi (3.24) 

with /3 G Q a rational number. The superpotential determined by (|3.24|) is then 



W{(j)) = A{(j) + i(3) . (3.25) 

One can then verify that under (p —(f) with the fiuxes inert, V{(j)) is invariant. In fact, 
this moduli-only transformation sends 

A 

W-^-—W. (3.26) 

Thus we find that it acts as a combination of an antiholomorphic transformation and 
a Kahler transformation ( |2.13|) , and hence when combined with spacetime parity it is a 
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symmetry of the low-energy theory. Note that when /3 = 1 , (|3.17|) holds and the additional 



Z2 (|3.19|) is also present. 

We can estimate the number of vacua with the complex conjugation Z2 symmetry by 
restricting the sum over at height H in ( p.lOj) to (p with vanishing real part. The number 
of such (p goes as H instead of H^, so the total number of vacua on the imaginary axis 
should go roughly as L log L. 

Let us now carry out a more careful analysis to check this rough asymptotic estimate. 
We are now looking at an overconstrained system, where extra conditions are imposed 
which geometrically constrain us to a hypersurface in the cone C over the space of fluxes. 
Thus, the continuous methods of [^| cannot be applied in this case and we must resort to 



a more discrete method of analysis. This is the first example we encounter in this paper 
where more number-theoretic methods are required. While the following analysis may 
seem somewhat particular to this problem, we will find parallel questions arising later in 
the paper, for example in the discussion of W = vacua on the torus in section §4.3. 

We can characterize the vacuum counting problem as a problem of counting integers 
subject to a certain set of constraints. We can write fluxes leading to Re (p = as 

A = k{p + iq), B = il{p + iq), (3.27) 

for k, l,p, q integers. To avoid overcounting, we can choose p + iq to be primitive in the 
sense that p, q are relatively prime. This gives a 1-1 correspondence between solutions of 
the vacuum equations with Re (p = and A^Aux < and integers kj,p,q satisfying the 
conditions that gcd(p, q) = I and 

kl{p^ + q^)<L. (3.28) 



We can write this as a sum 



iVvacua(i^;Re(^ = 0) = - Y. ^(^) (3-29) 

ip,q) = l n=l 

where d{n) is the number of divisors of n and counts the factors kl. The 1/2 arises 
because we want to fix S'L(2, Z) by taking (3 = l/k > 1. We can approximate the sum 
over p, q by an integral, where just as in (|2.29| ) we include a factor of ^(2) = 7r^/6 in the 

32 



denominator to remove integers with common factors, and where we use the use the fact 
that asymptoticaUy 

N 

^d{n) ^ N log N , (3.30) 

n=l 

which comes from the fact that asymptoticaUy the number of divisors of n goes as d{n) ~ 
log n. We then have 

iV,,,,,(L; Re = 0) = dp -r^^^og ^ ~ |-L(log . (3.31) 

Jl 4C(2) p^ Stt 

This gives an estimate for the number of vacua with the discrete Z2 complex conjugation 
symmetry. We have confirmed this result numerically by direct computation. Note that 
the naive asymptotic analysis above missed an extra factor of log L coming from the 
number-theoretic aspects of the problem. We have found that the fraction of vacua with 
this particular discrete symmetry goes as 

iVvacua(i^;Re(/) = 0) 9 (logL)2 



(3.32) 
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4. Flux vacua on = (T 

The torus is the simplest possible manifold on which to compactify string theory, so it 
is naturally one of the first places to study problems associated with flux compactiflcation 
We now turn to a study of vacua on a simple class of toroidal compactiflcations, 
namely those which can be described as a direct product (T^)"^ of 3 tori with equal modular 
parameter r. This example, while more complicated than the rigid CY, still has the 
advantage that the superpotential is of flnite order in the moduli. We solve the DW = 
equations for general fluxes on the symmetric torus, and show that the periods always 
take values in a fleld extension of degree at most 6 over the rationals Q. We study W = 
vacua and vacua with Z2, Z3, and complex conjugation symmetry. In our analysis we 
encounter examples of logarithmic L-scaling associated with a special cut through moduli 
space, nontrivial W = vacua, and interesting number-theoretic structures. We flnd in 
particular that the number ofW = vacua can be expressed in terms of class numbers, as 
anticipated by Moore , and that although vacua with discrete symmetries lie at special 



loci in the moduli space, they are reasonably plentiful. 

We parametrize the torus as follows (our conventions are as in 0). For completeness 
we will set up the system with a general T^, then specialize to the more restricted (and 
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hence more tractable) models of interest. Let a;*, for z = 1, 2, 3 be three pairs of periodic 
coordinates = + 1, j/* = + 1. Choose three holomorphic 1-forms dz'^ = dx'^ + 
r'^^dy^. The complex structure is completely specified by the period matrix r*-' . Choose 
the orientation 

dx'^ A dx'^ A dx^ A dy'^ A dy"^ Ady^ = 1 . (4.1) 



/ 



Then a convenient basis for H^{T^, Z) is given by 



Q!° = dx^ A dx"^ A dx^ , 
1 
2 



aij = -eumdx^ A dx"^ A dy^ , 



/5*' = -2W2/'Arfy"^A(ix\ 



(4.2) 



1 

— ( 
2 

/?° = dy^ A A dy^ . 
This basis satisfies the property 

j aiAp-^ = 6j. (4.3) 
The holomorphic three-form can be expressed in terms of the holomorphic one-forms as 

n = dz^ A dz^ A dz^ , (4.4) 

with the expansion 

n = a° + aijT'^ - (cof T)y + (3^ (detr) , (4.5) 

where 

(cofr),,- = (detT)r-i'^ = ^e,fc^e,p,T^V"^^ . (4.6) 

The modular group Q for a general flux compactification would include the SL{6, Z) 
group which acts on the toroidal moduli and the SL{2, Z) S-duality group. For the case of 
interest to us, where the is a product of three identical two-tori with complex structure 
parameter r, this will simplify to ^ = SL{2, Z)^ x SL{2, Z)^-. 

4.1. Superpotential and equations for flux vacua 

General fluxes can be expanded (with {27r)'^a' = 1) as 

(4.7) 

H3 = c°a° + c'^aij + dijP'^ + do/?° • 
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Here a^, . . .do are all integer parameters. The Gukov-Vafa-Witten superpotential can 
be expanded using ( ^T^l ) and ( [4.3[ ) as 



W={a^- (j)c'^) detr - {a'^ - 4>d^) (cofr)^^- - {hj - (j)dij) r'^ - {bo - <pdo) . (4.8) 

At this point, we restrict our discussion to the symmetric cases which will be of interest 
to us. The most general fluxes we will consider will take the form 

a'^ = a5'\ b,j = b5,j, c'^ = c5'^ , dij = d5,j . (4.9) 

That is, all flux matrices will be diagonal and proportional to the identity matrix. In 
addition, we allow general values of the integers aP ,bo,c^ ,do- 

The most general discrete solutions of the DW = equations with the fluxes ( [4.9| ) 
take the form r*-^ = r5*-^ - that is, the six-torus becomes a (T^)"^ where all two-tori have 
the same complex structure (as described for the W = Q case in f^). In fact, six-tori of 
this form have an enhanced symmetry 5*3 which permutes the three two-tori; the fluxes 
( [1.9| ) are the most general fluxes consistent with this enhanced symmetry. Moduli which 
would break the symmetry are by deflnition charged, and so as in ||^ we can argue that 
if we ignore the charged moduli and flnd solutions for r on the symmetric locus, these 
solutions will lift to solutions of the full theory including the charged moduli. 

For the symmetric torus (T^)'^ with fluxes ( [4.7|) , the superpotential (|4.8|) further sim- 
plifles to 

Vr((/>,r)=Pi(r)-(/>P2(r), (4.10) 

where -Pi(t) and P2{t) are cubic polynomials in the complex structure parameter r deflned 
by 

Pi(r) = a°r^ - Sar^ - 3br - bo , 

(4.11) 

P2(r) = c^T^ - Sct"^ - Mr - do . 



The Kahler potential for the theory truncated to r and 4> is given, using (|4.5|) and ( |2.2|) , 
by 

/C = -3 log(-z(r - r)) - log(-z((/> - ^)) , (4.12) 
and the D3-brane charge contribution for the symmetric torus is 

iVflux = ^oc*' - a°c^o + 3(6c - arf) . (4.13) 
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After some simplification, the vacuum equations take the form 

Pi(r)-^P2(r) =0, 

Pi(r) - 0P2(r) = lir-r) {P[{r) - <pP^{r)) . ^^'^^^ 

We can solve for the dilaton (p using one of these equations (when P2{t) 7^ 0), and then 
eliminate it from the other equation, leaving a single equation which is quintic in r, r: 

g(r, r) =P2(r) (Sar^ + 66r + 36o - r(3a°r2 - 6ar - 36)) 

(4.15) 

- Pi(r) (3cr2 + Qdr + 3cio - r(3c°r2 - 6cr - 3rf)) = . 

While rather complicated, this equation can be solved exactly. The complete solution is 
described in subsection §4.3. Before analyzing the complete system, however, we find it 
instructive to first consider a simpler model in which we only turn on a subset of the 
fluxes. In this simplified problem, which we analyze in subsection §4.2, the solution of the 
equations is more transparent, and already leads to interesting lessons about torus vacua. 

4-2. Vacua on special symmetric (T^)^ 

We can even further simplify the full symmetric torus problem by restricting the fluxes 
to have the special form 

d = —a , b = c = a^ = do = 0, bo.c^ arbitrary. (4.16) 

This reduced set of fluxes was discussed in §8.1 of 0, where restrictions were also placed 
on the allowed moduli. Here, we analyze all possible solutions to the equations of motion 
( [4.14| ) given this set of fluxes, with no restriction on moduli. 

The equations of motion in this case are fairly simple, but already lead to interesting 
phenomena including residual flat directions, nontrivial W = solutions, and an asymp- 
totic number of vacua with anomalous scaling. It will turn out that the set of all isolated 
vacua (those not part of a flat direction) coincides with the ones with restricted moduli 
considered in 0, and it is these that we will count. 

With the restriction (|4.16|) the tadpole condition ( [4.13|) becomes 



6oc° + 3a2 < L . (4.17) 



36 



4-2.1 Vacuum equations 

Consider the quintic ( |4.15| ). Imposing the flux restriction ( [4.16| ), and writing r = x+iy, 
we get two quintic equations in the real variables x,y. These equations have factors of 
y'^^y respectively. Since we are interested in solutions with y 7^ 0, we drop these factors 
and get the cubic and quartic equations 

2ac^ y"^ X = \—2a(Px'^ + 3a^ — hQ(P\ x = p2{x)x , 

(4.18) 

ac^ y'^ = ac^x^ + (^oc° ~ 3a^)a;^ + aho = Pa{x) . 

The first equation clearly has a solution when x = 0. In this case the second equation is 
satisfied when either a = or 

y' = 'i- (4.19) 

When a; 7^ 0, we can factor out an x from the first equation, and combine the two equations 
to get 

= p2{x)'^ — Aac^ p4^{x) 

(4.20) 

= 9a^ - IQaHoc^ + hl{c^f . 

Summarizing, we see that there are several continuous families of solutions, and one discrete 
family of solutions with the restricted class of fluxes we are considering. The continuous 
families of solutions arise when either a = 0, in which case we have an arbitrary imaginary 
value for r, or when ( |4.20|) is satisfled, in which case the real part of r is arbitrary, and 
the imaginary part is determined by 2ac^y^ = p2{x) (of course, if a = or c° = 0, r is 
completely unconstrained) . 

The appearance of these continuous families is not surprising. A certain number of 
fluxes must be turned on to stabilize all moduli, and in these cases we simply have too 
few nonzero fluxes to stabilize r completely. Similar flat directions for flux vacua have 
been discussed in e.g. mostly in situations where the fluxes preserve extended 

supersymmetry. 

The vacua we are particularly interested in are the isolated vacua, which we have 
found lie at the points 

r = i(^)"\ (4.21) 

where boc^ > 0. From (|4.14|) , we have 

, Pi(r) -3ar2-6, 



P2{r) c^T^-SaT 
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r. (4.22) 



Precisely these vacua were considered in §8.1 of H. 



Ji-.2.2 Counting vacua 



In enumerating vacua we have to make sure we are not double counting, so we must 
understand the action of S'L(2,Z)0 x SL{2,Z)r on the fluxes; the latter can be deter- 



mined using dz = dx + rdy along with the definitions (|4.2|) and ( ^77|) . The SL{2,Z)fj, 
transformation —1 can be fixed by assuming bo,c^ > 0. One finds that although more gen- 
eral transformations of each SL{2^Z) factor individually violate ( [4.16|) , the combination 
S-duality 

(j) r ^ -1/r (4.23) 

preserves it with the induced action on the fluxes 

6o ^ c° , a fixed . (4.24) 

One can deal with this redundancy by counting only vacua with > c^, which restricts 
( [4.21|) to the usual fundamental domain. 



To determine the scaling of the number of vacua in this family, we can use a continuous 
argument in which we ignore the discreteness of the fluxes, as the set of allowed fluxes 
forms a cone in the full 3-dimensional flux space we are considering. The boundaries of 
the cone are the conditions < c*^ < 605 and we want to count points in this cone with 
60c'' + 3a^ < L. In the continuous approximation, the number of allowed fluxes should be 

iVvacua(i^) ~ - / dhj dc' 2 J ^ , (4.25) 



1 



where the integration over a was just the integral of 1 from — a/(L — bod^) /3 to 



■\/(L — boc^)/3, and we have dealt with the redundancy ( 4.24|) in this approximation sim 



ply by dividing by 2. The bounds on bo, are determined by ( [4.17| ) and the requirement 
that ( [4.21| ) give only non-singular solutions, which rules out either vanishing. Performing 
the integration yields 

iVvacua(i^) ~ ^L^/^logL. (4.26) 

The logarithmic correction to power-law scaling for the number of flux vacua we find here 
is notable. Unlike the logarithmic term we found in §3.5.2, this logarithm is an example of 
the phenomenon mentioned in §2.3 in which the cone C containing the region of interest 
in moduli space intersects with the vanishing locus Z of the quadratic form A^Aux on the 
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sphere S'^^^~^, and goes to infinity. In this case, while the total volume of the cone for the 
full torus is finite, we have taken a slice through this cone of infinite area. The divergence 
of this infinite area is cut off by the lattice as / 1/x — > ^iLi 1/^ ~ log A^. Because this 
divergence comes from a region where the cone approaches the vanishing locus Z, in fact 
the continuous approximation is not really valid in this case and we should check the 
overall coefficient by performing a discrete calculation. It turns out, however, that in this 
case imposing the lower cutoff on the fiuxes at 1 precisely gives the right coefficient, so the 
answer above is correct. 

4.2.3 W = vacua 

For the vacua ( |4.21[ ) , ( |4.22| ) the superpotential takes the form 



W{<P = T = i{bo/co)'/') = 2\^i3a - ■ (4.27) 

V c" 

W = precisely if Bqc^ = 9a^. This means A^Aux = 120^. The number of W = vacua 
with Afflux < is then given byS 

A^vacua(i^;Vl/ = 0)= Yl c^(9a2)~— ^-(log(L))2 + 0(v^logL). (4.28) 

a=l 

Therefore, we see that the fraction ofW = vacua in this simple toy model is 

iVvacua(i^; W = 0) 9 logL 

iVvacua(i^) ^ 87r2 L 



(4.29) 



4.2.4 Discrete symmetries 

Two enhanced symmetry points of ^ = SL{2, Z)(j, x SL{2, Z)j. could in principle arise 
for solutions with r = (p = i and t = (p = exp(27rz/3), but only the former is consistent 
with ( [4.21| ); this is the fixed point of the Z2 symmetry action discussed previously, arising 
for fiuxes bo = c^. 

As before, we can ask whether this transformation of the moduli and fiuxes descends 
to a symmetry only acting on the moduli when the vacuum sits at the fixed point. Since 
the combined transformation is always a symmetry, this is equivalent to asking whether 



^ Thanks to N. Elkies for explaining to us how to derive the constant factor here. 
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transforming the fluxes only is a symmetry, as in (|2.37| ). It is obvious from ( [4.24 ) that at 
bo = c° this is indeed the case. Checking exphcitly, one flnds that 

W{ct>, t)^W{-^,--) = -^W^(<^, r) . (4.30) 



Checking the action on the Kahler potential ( [4.12|) as well, we flnd that this transformation 



is indeed a Kahler transformation on the low-energy theory, and hence a symmetry. 

To count the enhanced symmetry points, we should count the integer choices of bo, a 
satisfying 

bl + 3a^ <L . (4.31) 
The continuous approximation is valid here 



-^vacua 



(L,r = <^ = ^)~2^ dbo ^J{L-bl)/?,^^L . (4.32) 



The ratio of vacua with enhanced discrete symmetries to generic vacua therefore scales like 

iVvacua(-^, T = (/) = i) 3% 1 



A^vacua(i^) 4 yZlogL ■ 

4-3. Symmetric (T^)"^ 



(4.33) 



Now we consider the more general symmetric torus vacua characterized by vacuum 
equations ([4.14| ) , ( [4.15| ) with arbitrary flux parameters a*^ , . . . , (io ■ As for the special sym- 
metric torus with restricted fluxes considered in the previous subsection, we can again solve 
the equations completely, although the details are somewhat more complicated. We have 
used the exact solution to numerically generate a large number of sample vacua, which 



we can compare with a prediction of Denef and Douglas [^. We can also give a simple 
algebraic characterization of = vacua in the general case, which allows us to compute 
the asymptotic fraction of flux vacua in this category. We also discuss discrete symmetry 
points in this model. 

4-3.1 Vacuum equations 

Using again the parameterization t = x + iy, we can write the quintic ( ^4.15| ) as a pair 
of coupled equations in x, y of the form 

gi(x)y2 = qsix) , 

(4.34) 

qo{x)y^ = q4^{x) , 
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where qi{x) are respectively quartic, cubic, linear, and constant functions of x for i = 
4,3,1,0 with coefficients which are polynomials in the flux parameters a*^, ...,(io with 
integer coefficients. These functions are given explicitly in the Appendix. Combining 
these equations, we have 

qi{xfq4{x) = qo{x)q^{xf . (4.35) 

Naively this looks like it should be a 6th order equation in x, but remarkably the terms of 
order 4, 5, and 6 all cancel, giving us a cubic equation for x, 

asx^ + a2X^ + a\x + a'' = . (4.36) 

The coefficients in this cubic equation are presented explicitly in the Appendix. Cubic 
polynomials can be solved explicitly in closed form, so for any set of fluxes we can compute 
the values of x satisfying ( [4.36| ). When either q\{x) 7^ or q^ix) 7^ we can then solve for 
y. As we found in the previous subsection, certain flux combinations do not fix the modular 
parameter r uniquely. For example, ([4.36| ) vanishes identically for the flux combination 



a" = 0, a = 8, 6 = 4, 60 = -2, c° = -1, c = 7, = 3, rfo = -2 , (4.37) 

and thus x is undetermined for these fluxes, indicating the presence of a flat direction 
in the complex structure moduli space. We are primarily interested in vacua where all 
complex structure moduli and the axio-dilaton are fixed, so we do not bother trying to 
classify the fluxes with this property. 

In general, we see that x belongs to a field extension Tx of degree at most 3 over 
Q. When q\{x) 7^ 0, y obeys a quadratic equation with coefficients in T^- The resulting 
T = x + iy lives in a field extension J-' of degree at most 6 over Q. For any particular values 
of the flux parameters a*^, . . . , (io we can solve the equations explicitly and give algebraic 
expressions for the possible allowed values of r. 

4.3.2 Counting vacua 

In Denef and Douglas used the approach of to estimate the asymptotic 



number of vacua with A^Aux < L on the symmetric torus as L ^ 00, including a proper 
treatment of the 5'L(2,Z)^ x SL{2, Z)j. gauge fixing. They found that the total number 
of vacua should go as L^, and that the distribution on the complex structure modulus 
r should be proportional to l/(Im r)^. We have used the exact solution described in 
the previous subsection to perform an independent test of the number and distribution 
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of vacua on the symmetric torus. We used a Monte Carlo algorithm to generate a large 
number of solutions, gauge fixing SL{2, Z)^ with conditions on the fluxes, and gauge fixing 
SL{2^Z)r by only keeping solutions with r in the fundamental domain. Our numerical 
results are compatible with the results of , both for the total number of vacua and the 
distribution on moduli space. We found the L'^ scaling for L in the range 20 < L < 100, 
where we have A^vacua ^ CL'^ with C ~ 0.13. We have also confirmed numerically that 
when Afflux > 20, the distribution of vacua is compatible with l/(Im r)^. For values of 
L < 20, the number of vacua shrinks more rapidly than as L decreases, indicating that 
below this value the number of vacua is suppressed due to lattice effects. At small values 
of L the existing vacua are dominated by small values of Im r, indicating that the region 
of moduli space with large Im r is near a void in the lattice structure. 




We have graphed the distribution of vacua in the r fundamental region for small values 



of Afflux ^ 18 in fig. 2. As in the rigid CY model, for which a graph appears in [26|, there 
is quite a bit of structure to the points. In particular, something like the void found in 
||26|| appears in the vicinity of the imaginary axis. This phenomenon may be related to 
the large number of vacua which arise along this axis, including those described in §4.2, 
many of which have a discrete symmetry under r — > — r, as we discuss in §4.3.4. This void 
structure may also be related to the fact that the imaginary axis is one place admitting 
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one-parameter families of fliix vacua, as discussed in §4.2. Since we are now considering 
vacua with periods which generally are associated with a field extension over Q of degree 
up to 6, the allowed values of r have much greater freedom. It would be nice to develop a 
better understanding of the structure of this family of vacua. 

Some comments on the details of our numerical investigation may be of help in future 
efforts of this type. We gauge-fixed the SL{2, Z)j. symmetry by choosing a'^ = 0, < c < a. 
This gauge-fixing is not complete, as it does not deal with situations where c'^ = or a = 0. 
We did some further gauge fixing to deal with most of the residual gauge symmetry, but 
these cases should not contribute at leading order in L. For each value of A^Aux < Lwe then 
chose random fiuxes, selecting 5 of the nonzero fiuxes in a range —L < d, . . . < L, and 
determining the remaining two (c°, bo) by choosing all possible factorizations of the value 
6oc° compatible with the other fiuxes and ( [4.13| ). (This step increases the effectiveness of 
the algorithm by a factor of since we essentially test values of the fiuxes at once.) 

It is worth commenting on the range in which we tested the fiuxes. While a gauge 
fixing of SL{2, Z)^ by choosing ^ in the fundamental domain should give fiuxes of order 
-\/L, gauge fixing in fiux space generically gives fiuxes of order L. This is easy to check in 
the rigid CY model, where the fiuxes are uniformly distributed in a range < H < L. For 
the torus model we are considering here, empirical tests at a range of values of L, including 
a systematic check for small L < 9 and a Monte Carlo check for L < 100 found no vacua 
with individual fiuxes > A^Aux, while vacua do arise with individual fiuxes which saturate 
the upper bound. 

4.3.3 W = vacua 

Now let us consider the problem of finding vacua with W = on the symmetric torus. 
It follows from ( [4.10| ), ( [4.11|) , and ( [4.14| ) that to find W = vacua, the modular parameter 
r of the tori must solve the two cubic equations 



(4.38) 



P2(r) = c°r^ - 3cT^ - Sdr - do = . 



(4.39) 



Prom the second equation of ( [4.14]) , we then have the dilaton 




(4.40) 



The contribution to the D3-brane charge is again (|4.13| ). 



43 



In order to be physical, r must have a nonzero imaginary part. But since the coef- 
ficients of (|4.38| ) and ( |4.39| ) are real (and in fact integral), complex roots must come in 
conjugate pairs. Therefore, ( |4.38| ) and ( [4.39| ) must share a common quadratic factor. In 
fact as argued in 0, they must factorize over the integers so that 

Pi{r) = {fT + g)P{T), (4.41) 

P2{r) = kP{T), (4.42) 

where 

P(r) = Zr^ + mr + 71 , (4.43) 

with f , g.k.l.m.n G Z. Here we have used the SL{2,Z)ff, symmetry to set c° = and 
< —^0 < ~dQ, making P2{t) = a quadratic equation. 

This system of equations is closely parallel to that encountered in §3.5.2. There we 
had a pair of fluxes in Z[z] with a common factor. We can use a very similar analysis 
to describe the number of = vacua in this case. Before going through the detailed 
analysis, however, we can use the method of heights to give an estimate of the scaling 
we expect. In this case, the height is simply the scale of the integer coefficients in the 
various polynomials. As mentioned in §2.4, the height of a product of polynomials goes 
as the product of the heights of the factors (this statement is transparent for polynomials 
with integer coefficients). We wish to know what fraction of the pairs of cubic polynomials 
( [4.38|) , ( [4.39|) share a common quadratic factor. The number of quadratics with coefficients 
of height q goes as . We multiply this quadratic by two linear functions, each of height 
< y/Ti/q, as the coefficients in the cubics should be of height \fL. The number of such 
linear functions is L/q^. Thus, we expect the total fraction of cubics with coefficients of 
order \/L which admit a common quadratic factor to go as 

Furthermore, we expect the total contribution from quadratics of height g to go as 

Note that while we have used the height here as a property of the common quadratic 
polynomial factor P, there is a close connection between this analysis and the algebraic 
properties of the modular parameter r where P{t) = 0. The set of r which satisfy quadratic 
equations with integer coefficients are the quadratic imaginary integers, which lie in com- 
plex extensions of degree 2 over Q. By summing over all polynomials of all heights, we are 
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equivalently performing a sum over all complex degree 2 extensions of Q and over all r 
which lie in each of these extensions. As we see in the following analysis, the contributions 
to the sum naturally group themselves in terms of the field extension, which is determined 
by the discriminant of the quadratic -P(t). 

Now let us go through the analysis in more detail, following the same approach as 
in §3.5.2. As there, we fix the redundancy in the parameterization ( [4.41|) , (|4.42|) , ( [4.43| ) 
by choosing gcd(/, m, n) = 1; -P(r) is thus a primitive quadratic over the integers. We 
furthermore fix the sign ambiguity by setting / > 0. From the coefficients in ( [4.38| ) and 
( [4.39| ), we see that there are further mod 3 constraints on the integers, so that 

kl = km = gl + fm = gm + fn = (mod 3) . (4.45) 

In terms of the parameters f,...,n, the tadpole condition becomes 

< {fk){4ln - m^) = 3Nfi^^ < 3L . (4.46) 

We still need to fix the redundancy associated with SL{2, Z\. In this case, it is 
convenient to do this by simply solving P(r) = and requiring that one of the roots lie in 
the fundamental domain. From 



— m ± Aln — m^ 
r = ^^ . (4.47) 

we see that the necessary conditions are 

-l<m<l<n, (4.48) 

with m < when / = n. These conditions guarantee that the discriminant is negative, or 
that 

= 4/n - > . (4.49) 
In addition, the previous 5'L(2, Z)<^ fixing requires 

0<g <k . (4.50) 

Note that the SL{2, Z)t- fixed points of r = i, r = exp(7rz/3) correspond to fiuxes living in 
an SL{2, Z)rj, conjugacy class which is invariant under the action of the relevant SL{2, Z\ 
generators. For example, consider P\ = r"^ + r, P2 = + 1. Under the action oi Sr '■ t ^ 
we have Pi —P2,P2 -Pi, which can be put back in canonical SL{2, Z^^ form 
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by S(f,, which returns us to the original polynomials. The consequence of this is that we 
weight all vacua, including those at the Z2 and Z3 points, with weight 1. 

To compute the total number of vacua at fixed A^Aux, we now simply sum over all 
combinations of integers satisfying r G J?-d (the fundamental domain), ( ^.45| ), ( |4.46D , and 



( [4.50|) . Aside from the mod 3 constraints, this is a very similar problem to that encountered 
in §3.5.2. Let us first consider a case where the mod 3 constraints are simple to incorporate, 
namely when il>(/,m, n) 7^ (mod 3). In this case, the constraints (|4.45|) simply force the 
integers /, g, k to all be multiples of 3. We write / = 3/, g = 3'g, k — 3k. We now can write 
the sum for these values of D as (using r = fk) 

l,Tn,n:TeJ^,S\D r=l -^|^ 

The factor of k inside the summation includes all possible values oi g < k (including 
(7 = 0). This summation can be written in terms of standard number-theoretic functions. 
The number of relatively prime integers /, m, n giving rise to a r in the fundamental domain 
with fixed D is known as the class number h{D). For a discussion of class numbers in the 
physics literature in a related context, see [^. The sum over k is just the sum over divisors 
of n. Thus, our summation in the case where 3\D reduces to the number-theoretic sum 

J2 h{D) J2 • (4-52) 

DeSZ n=l 

The mod 3 relations (f4.45| ) only make a small change to this analysis. Consider for example 
the case where 3\D because I = m = (mod 3). In this case the mod 3 constraints only 
force / = 3/, while g, k can be arbitrary. The summation above is now exactly the same 
except that instead of 3D in the denominator of the upper bound for n, we just have D. 
One can check that in all cases where D = (mod 3), the correct summation always has 
the factor of D in place of 3-D. Thus, we can write our final formula for the number of 
W — vacua on the torus 

DeZ n=l 



where 



m = {o;. if 15^ (mod 3); 
' 3U, otherwise. 
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The fact that the total number of vacua in this case can be expressed in terms of the 
class number and other number-theoretic quantities was anticipated in [|5^, and seems 
indicative of a deep relationship between W = vacua and arithmetic. 

We can furthermore separate out the contributions to the total number of W = 
vacua ( [4.53|) based on the value of the modulus r, given by the solution (|4.47| ) to the 
primitive quadratic parameterized by relatively prime integers (/, m, n): 



t{Aln — m^) J- 2 

' It 



k 

For example, the number of vacua at r = i, associated with the quadratic -P(r) = + 1 
and the integers m, n) = (1, 0, 1) is 7r^/1728 ~ 0.00571 L^, while the number of vacua 
at r = exp(7rz/3) {P = - t + 1) goes as ttVIOS ^ 0.09139 L^. 

From ( |4.55| ) we see that, just as calculated using the height argument, the contribution 



from each primitive quadratic goes asymptotically as L^, with a coefficient which is of order 
l/g^. The primitive quadratics at height q combine to give a total contribution at height 
q of l/q'^. Using ^ 1/n^ = 7r^/6, a rough estimate for the total number of vacua is then 
given by 0.097 (tv^/G) ^ 0.159L2. 

To determine the precise asymptotic coefficient of in the total number of vacua, 
we now simply need to carry out the sum 

This sum converges quite rapidly. The contributions at fixed q = n seem to go as 
indicating that more subtle number theoretic considerations correct the heuristic argument 
above by log(log q) terms. Note also that the contributions at different q jump about with 
a periodicity of 6 due to the divisibility properties of the integers. 

In fig. 3 we show the exact number of = vacua with A^Aux < L and the curve 
representing the asymptotic form at large L in a log-log plot. Note that while the sum 
( [1.56| ) converges rapidly in q, the exact value of A^Aux approaches the asymptotic form 



much more slowly. 
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Fig. 4: Distribution oi W = vacua in complex structure fundamental domain 
for L = 5000 for small values of Im r. 



In fig. 4 and fig. 5 we illustrate the distribution of = vacua in the complex 
structure fundamental domain for small and large values of Im r respectively. The curves 
visible in fig. 4 can be interpreted as the loci of values of r associated to vacua with a linear 
constraint imposed on m, n); using ( |4.47|) , one sees that constraints not involving n will 
fix X = Re r, while for n = Cim + C2I one relates x to y = Im r as + = —2Cix + C2. 

The large scale peak structure in fig. 5 can be deduced from ( [4.47|) with the additional 
modular group fixing constraint (|4.48| ). In particular, the heights of the maximal peaks 
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Fig. 5: Distribution of VF = vacua in complex structure fundamental domain 
for L = 2000 for large values of Im r. 



are given by 

max(Im r) ~ — at Re r = 0, ±0.5 . (4.57) 

The next peaks have height Im r ~ V^/4 at Re r = ±0.25. We also confirm numericaUy 
that the distribution of vacua in the complex structure fundamental domain is in accord 
with l/(Im r)^. 
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Fig. 6: Void structure of distribution of Vl^ = vacua in dilaton fundamental 
domain for L = 600. 
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In fig. 6 we illustrate the distibution of W = vacua in the axio-dilaton fundamental 
domain, where the dilaton following ( [1.40| ) is given by 

. fr + 9 



k 



(4.58) 



which we map to the fundamental domain. The void structure that appears in this case 
has an analogous explanation to the voids which appeared in the rigid CY case. 

4.3.4 Discrete symmetries 

Two enhanced symmetry points arise for (T^)^, related to fixed points of the modular 
group Q = SL{2,Z)fj, x SL{2,Z)t-: (j) = t = i and (p = t = exp(7rz/3), with diagonal 
SL{2, Z) transformations of the dilaton and complex structure. We also have a complex 
conjugation symmetry which takes r — r, (p —>■ —(f) similar to the symmetry of the rigid 
CY discussed in section §3.5.2. We discuss each of these discrete symmetries in turn. 

Z2 C SL{2,Z) 

For the point (p = t = i there is a combined Z2 discrete symmetry 

r-^r' = --, (f) (/,' = --. (4.59) 
r (p 

This symmetry arises for both W = and generic no-scale fiux solutions at the fixed point. 
From ( ^.14| ) we determine the fiuxes leading to vacua at this fixed point: 



a° = -do, 6o = c°, a = -d, b = c. (4.60) 
The superpotential determined by fiuxes ( [4.60| ) then transforms under ( [4.59| ) as: 

W{t, c/>) ^ W{r', (/>') = 47^(^' (4-61) 

We can count the number of vacua at (p = t = i easily by simply imposing the corre- 
sponding linear conditions on the fiux. Since the moduli have heights of order 1, we expect 
the number of vacua to scale as with a coefficient of order 1 when we only impose the 
DW = conditions (|4.14| ), and as L when we also impose the W = condition. We have 

Afflux = + 36^ + (a0)2 + bl < L , (4.62) 

so the total number of vacua is 

iVvacua(i^; <P = T = i)r^ —L^ . (4.63) 
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Imposing the further constraint W — gives us the further relation on the fluxes 

a" = -36, 60 = 3a , (4.64) 

leading to 

Afflux = 12(a2 + 62) < L . (4.65) 
Then the total number of vacua with 1^ = and the discrete Z2 symmetry is 

iVvacua(i^; (j) = r = i;W = 0)^^L. (4.66) 

ZsCSL{2,Z) 

For the point 4> — r — exp(7rz/3) there is a combined Z3 R-symmetry 



1 , „ 1 

r — r : 



r-1' ^ ^ (h-1 



(4.67) 



In this case the symmetry arises only for W — solutions at the fixed point; the fiuxes 
leading to these vacua are 

a^^-2b + d, bo^c^^-b-d, do ^ b - 2d , a ^ -b , c^-d. (4.68) 

The superpotential for such fiuxes as a function of the complex structure and dilaton 
transforms as 

W{t, <P) ^ WiT\ <P') = W(r, <P) . (4.69) 

In particular, at the fixed point W — exp(7rz/3)VF which requires W = 0. We expect 
again that the number of vacua at this point will go as L with a constant coefficient of 
order 1. The value of A^Aux given the fiuxes ( [4.68| ) is 

Afflux = 3(6^ -bd + d^) < L. (4.70) 

Hence the total number of vacua with = and the discrete Z3 symmetry is 

iVvacua(i^; </> = r = exp(7rz/3); = 0) ~ . (4.71) 



51 



Complex conjugation Z2 

The torus possesses a complex conjugation symmetry C of the form ( p.38|) , with U = 
diag(— 1, 1, —1, 1); the fuU action is 

a*^ — > — , h ^ —h , c—^—c, do —do , r ^ —r , (p —(p . (4.72) 

For certain combinations of fluxes, the action of C takes the fluxes to themselves, up to a 
possible overall sign. For such fluxes, the transformation r — r, (p —cp acting on the 
moduli alone describes a symmetry of the low-energy theory when combined with four- 
dimensional parity; as for the complex conjugation symmetry discussed in section §3.5.2, 
this is an antiholomorphic transformation combined with a possible Kahler transformation 
(if there is a sign). 

The points in moduli space which are invariant under the action of this symmetry are 
those for which both cp and r are pure imaginary, so we can write 

(P = i\(Pl T = i\T\. (4.73) 

At such points, the equations (|4.14| ) take the simplified form 



rfo|</>| +a°|r|^ = 0, h = c\(P\\T\, 60 = c°|(^| |r|^ + a|r| = . (4.74) 

Prom (|4.34| ), where in this case x = 0, y = |r|, we see that is always rational, and that 
sometimes is rational, so that all solutions of ([4.74|) are associated with field extensions 
of degree either 2 or 4. 

Let us first take the case where the field extension is of degree 4, so that r = i{p/q)^^'^, 
with gcd(p, q) = I and such that p, q are not both perfect squares. In this case, |r| and \t\^ 
are independent over Q. For nonzero cp this means that the equations ( [4.74| ) cannot all be 
satisfied nontrivially: either all terms in the first two equations must vanish or all terms in 
the second two equations must vanish. This means we must set either a^ = b = c = do = 
or a = bo = c'^ = d = 0, with the nonvanishing fiuxes obey two relations with (p, r. 
The nonzero fiuxes in this case are precisely those which are taken to themselves up to a 
possible sign under the action of C. Thus, the vacua where (p, r are both pure imaginary 
and live in the field extension of degree 4 over Q are precisely those where we expect a 
low-energy symmetry. 

We can use a continuous fiux argument to enumerate these vacua, as the modular 
parameters are determined from the fiuxes, so we expect order vacua in this case. In 
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particular, in the situation where the aU terms in the first two equations vanish, we have 
a*^ = b = c = do = 0, leaving us with the fluxes a, 605 d satisfying the last two equations 
in ( [4 .741) . This is almost identical to the vacuum counting problem from §4.2.2, except 
that now we no longer impose the condition d = —a. This gives us an extra factor of a/L 
above ( [4.26 ). As in that case, however, we expect additional log L terms here also. 

Note that the structure of the constraint equations has simplifled dramatically in this 
set of vacua. In general, for an extension of degree 4, we would expect 8 equations for the 
8 fluxes. Because of our choice of fluxes and the structure of the equations which results 
when we choose moduli of the particular form i{p/q)^/'^, these equations have collapsed to 
only two constraints relating the remaining fluxes and the modular parameters. Applying 
the height-based arguments of §2.4 in this situation to, for example, the last two equations 
in ( [4.74|) , we see that the last equation suggests that if r is deflned by a quartic equation 
qr'^ + p = of height H = max (1^1, then (p/r has height V^, so the third equation 
gives a single constraint reducing the flux combinations by a factor of \/L ■ \/LH. There 
are order H r's at height H, each of which then contributes ^ L/H vacua, for an expected 
total of vacua distributed fairly uniformly across the height of r. 

In the case where the fleld extension is of degree 2, the situation is more complicated. 
An argument based on heights suggests that there are order such vacua. These vacua, 
however, do not have the complex conjugation symmetry of the low-energy theory. 



5. Vacua in simple Calabi-Yau hypersurfaces 

We turn now to the last class of examples, the set of one-parameter Calabi-Yau mani- 
folds defined as hypersurfaces in weighted projective space. We analyze these vacua at and 
around a particular point in moduli space, the Landau- Ginzburg point. We find that the 
structure of the number field in which the periods there take values plays a fundamental 
role in determining the nature of the solutions. We focus on the two models in which vacua 
can arise at the LG point, analyzing in each case the total number of vacua, the number of 
vacua with W = 0, and the vacua with discrete symmetries. We find that in one case the 
generic arguments in §2.4 are violated due to a convergence in the constraint equations, so 
that more vacua arise than naively expected. We find in general that Landau-Ginzburg is 
a very special point in the moduli space, where the number of vacua with W = and with 
discrete symmetries can be of the same order as the total number of vacua, when nonzero. 
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The Calabi-Yau manifolds of interest, for k = 5,6, 8, 10, are defined by the equa- 
tions 

4 

k = 5 : — bi/j xoa;iX2X3a;4 = , Xj G , 

i=0 
4 

k = 6 : + 2xq — Gi/j X0X1X2X3X4 = , x^ G WP'*i^i^i^i^2 , 

i=l 
4 

= 8 : xf + 4xq — 8tp X0X1X2X3X4 = , Xi E WP*i^i^i^i^4 , 

1=1 
4 

k = 10: ^x-° + 2xi + 5xo - 10-0x0X1X2X3X4 = 0, Xj G WP^i, 1,1,2,5 • 



(5.1) 



1=2 

The geometry of the spaces was studied in foUowing the seminal work of Candelas 
et al. on the quintic (see also ||61|,^ for discussions of the one-parameter models), 
and we follow their notation. The mirrors Wk of these spaces have one complex structure 
modulus, parametrized by if^ in (|5.1|), and are obtained by dividing out by the appropriate 
discrete group Qk [ |63| . In the spirit of 0, however, we work with the rather than 
their mirrors. As long as we turn on only ^fc-invariant fluxes, those complex structure 
moduli that are not invariant under Qk can be set to zero, as they are in ( |5.1| ), leaving us 
with effectively one-parameter models. This simplification is analogous to that made in 
the previous section to reduce to a model with a single r modulus. 

The moduli space of t/j for each is quite similar: all four have (mirror) Landau- 
Ginzburg points, conifold points, and large-complex structure points at -i/j = 0, 1, 00 re- 
spectively. For these one-parameter models, we focus on physics at and near the Landau- 
Ginzburg (LG) point -0 = 0, which is a fixed point of the modular group. Because we 
work with a perturbative expansion of the periods around this point, while we can solve 
perturbatively for solutions of the vacuum equations DW = 0, we cannot identify points 
in the moduli space with = or with discrete symmetries except precisely at the LG 
point, since these characteristics cannot be definitively identified in perturbation theory. 

The period vector of each of the at 'i/' = is associated with a different field 
extension over Q, and this algebraic structure is vital to understanding the kinds of vacua 
that exist (or fail to exist) there. For each of the models, we ask whether vacua exist 
at the LG point, whether W = vacua exist there, under what circumstances discrete 
symmetries arise in the low-energy theory at this point, and how many vacua with each of 
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these properties there are. Despite the naive similarity of the models, we shall see that Mq 
and Mg give rich sets of solutions on the LG locus, with accompanying discrete symmetries 
and = vacua in some cases, while M5 and Mio do not. In §5.4, we also count vacua 
in a small neighborhood around the LG point and the conifold point in Mg, to check the 
predicted V'^ scaling. 

Near the Landau- Ginzburg point '0 = the periods admit an expansion 



1 



expi 



{k-iyi 



< 1. Here a is the k^^ root of unity a = exp(27rz//c), 7 



(5.2) 



valid for 

and i^a, G = 0, 1, . . . 4 are the indices characterizing the weighted projective space. We have 
rescaled the holomorphic three- form O by 1/?/' relative to [|5^ . 

A convenient choice for the period vector is w'^ = {w2{i{^),wi{'il'),wo{^{j),Wk-i{i^))- 
The periods in a symplectic basis = {Qi, Q2, z^, z^) can be expressed in terms of this 
Picard-Fuchs basis by means of the linear transformation 11 = m/j ■ w where the matrices 



ijik have rational elements and are given in |5^. The symplectic basis period vector H 
then has the generic expansion 
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+ (9(^2) = Co + ci^ Pi + 0(^') , (5.3) 



where the Ci are complex constants and the pi are four- vectors with components that are 
rational linear combinations of powers of a. 



5.1. Vacua at the Landau- Ginzburg point 

We turn on flux vectors / = (/i, /2, /s, f^) and h = {hi, /i2, ^3, ^4)- The F-flatness 
conditions (^]§) at ■i/' = are then 



(/-(/>/i)-po = 0, {f-(j)h)-pi = 0. 



(5.4) 



In the cases k = 8 and /c = 10 a subtlety arises: as is visible in (|5.2|) , ci = in this 
case, and naively the D^W equation is trivial. However, another consequence is that the 
inverse metric g"^"^ diverges, such that the vanishing of g'^'^ D^W D:^W is nontrivial. The 
equations (|5.4|) should then be imposed with the vector p2 substituted for pi . 
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5.1.1 Vacuum counting and cyclotomic fields 

To solve the equations (|5.4|) , we must consider the algebraic properties of the vectors 
Pi. The Pi take values in rational combinations of powers of the k^^ root of unity a, which 
defines the cyclotomic field J-'{Q,)k- As we discussed in §2.4, the degree of J^k is given 
by the Euler totient function = (p{k) (|2.18| ). As is by now familiar, each equation in 
( ^■4|) becomes Vk constraints on the integer fluxes when the moduli are specified. Since 
there are 8 flux integers, we need 2(p{k) < 8 or (p{k) < 4 if we are to generically expect 
nontrivial solutions; for (f){k) =4 there are as many equations as fluxes, but one would 
expect vanishing fluxes to provide the unique solution. This is just the application of the 
logic that went into (p.20|) to this particular class of examples. 

For the various one-parameter models, one finds ^(6) = 2, but ^(5) = ^(8) = ^(10) = 
4. Hence according to the argument above, we expect nontrivial flux vacua at the Landau- 
Ginzburg point on Mg, but not for the other spaces. As we shall see in the next section, 
there are indeed LG vacua in the Mq model. However, we shall find that for certain 
restricted choices of the dilaton, there are LG vacua for Mg as well. 

These LG vacua at Mg exist despite the counting of equations because although ( |5^ ) 
produces 2(j){8) = 8 conditions on the fluxes, they are not all independent; for certain 
choices of the axio-dilaton (p, the equations degenerate into fewer conditions on the fluxes. 
This was one of the loopholes we anticipated in §2.4, that can invalidate the expected 
scaling ( |2.2C1| ) and yield larger numbers of vacua. Another example of this reduction in 
constraints appeared for the complex conjugation symmetry on the torus discussed in 
54.3.4. 



Returning to the logic of the analysis that led to (|2.20|) , we can try to parametrize our 
expectations for this nongeneric situation. If for a given choice of the axio-dilaton VEV 
there are m relations between the ry = 2^(62,1 + 1) equations on the integer flux parameters, 
as we shall find in the k = 8 case, then one should modify the definition of k in ( |2.20|) : one 
finds an effective Kef / = 263 — {rj — m), and the expected scaling of the number of solutions 
with L is A^vacua ~ L^'^^i'i'^^'^ . For the k = 8 model we shall find, for special values of 
that m = 4. For M5 and Mio, this phenomenon does not occur and there are no vacua at 
the LG point, as expected. 

The type of analysis we carry out in this section should be applicable for any Calabi- 
Yau which admits fiux vacua where the periods have ratios which are valued in an extension 
of the rationals of finite degree. While as we saw, for the torus all fiux vacua have this 
property, it may sound like an unusual condition for fiux vacua on a Calabi-Yau. But in 
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fact it has been shown to occur at the LG locus in a wide variety of multiparameter models 



|6J]. For example in one of the two-parametric CY models there is a flat direction ofW = 
vacua constructed in 0, with arithmetic periods — though in this case the dilaton is not 
constrained due to the flat direction. It may be that this sort of arithmetic structure for 
the periods is in fact a reasonably common occurrence for flux vacua. 

5.1.2 The modular group 

Before studying the one-parameter models, we must also consider the issue of gauge 
redundancy. The monodromy group F of the complex structure moduli space has two 
generators: A, which generates phase rotations "0 —>■ aij) with a = exp{2ni/k) around the 
LG point at "0 = 0, and T which corresponds to the logarithmic monodromy Q2 G2 + 
around the conifold singularity '0 = 1. A alone generates a C F subgroup, with an 
associated fixed point at '0 = 0; T on the other hand is of infinite order. 

In principle both SL(2, Z) and the modular group F will generate an infinite number 
of copies of any fixed vacuum. The former can be dealt with as it was for the previous 
examples, with a constraint either on the fiuxes or by choosing a fundamental region for 
the dilaton (p. We can try to fix F similarly. What is a fundamental region for the action 
of A and T? 

A, which rotates around "0 = by a k*'^ root of unity, obviously has as a fundamen- 
tal domain a wedge of 2n/k radians around this point. For T, on the other hand, by 
choosing the periods (^]^) we have already implicitly chosen a fundamental domain. T- 
transformations shift one between an infinite number of sheets defined around the conifold 



point '0 = 1, where the periods have a logarithmic singularity. By using the expansion ( |5.2| ) 
we have confined ourselves to a definite sheet. Hence T-images of the vacua we find will 
be solutions to the F-flatness equations with different periods, and we will not encounter 
them with our gauge choice (|5.2| ). 

Consequently, looking for vacua at or near '0 = using ( |5.2| ) we find only the k copies 
of a given vacuum generated by A; the T redundancy is already implicitly fixed. Numerical 
calculations with only SL{2, Z) fixed confirm this by finding only a finite number of vacua 
at fixed Afflux at '0 = 0, rather than the infinite number that would result if all T-images 
appeared. In the examples we study, the effective modular group is thus SL{2, Z) x Ti^. 

We now turn to a study of the one-parameter models in detail. Since only Mq and 
Mg have vacua at the LG point, we focus on those two cases, enumerating generic LG 
vacua as well as those with discrete symmetries and = 0. In the case of Mg, we also 
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report on a study of the number of generic vacua in a small region around the LG point, 
where solutions are found by using the expansion (|5.2|) and solving the equations using 
perturbation theory in t/j. 



5.2. Landau- Ginzhurg vacua for Mq 



Let us consider the case of /c = 6 first, the hypersurface in WP^J 1112- The root of 
unity is a = exp(7rz/3) = (1 + z-\/3)/2. Because = a — 1, the cyclotomic field jFg is just 
of degree two, with basis elements {1, a}. The vectors pi can be written 
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(5.5) 



Solutions of (|5.4| ) then exist as long as (p is chosen in J^q: 

(j) = ti + at2 , 



(5.6) 



with ti, t2 rational. Because the field extension in this case is of degree 2, the equations 
( ^.4|) for a fixed value of (p impose 4 real relations on the fiuxes ft, hi. These relations can 
be conveniently expressed in terms of the monodromy matrix A generating rotations by a 
root of unity around = 0: 

AIi{i;) = aU{atp) , (5.7) 

given explicitlyli by 

/I -10 1 ^ 

A = 



1 
3 -3 1 
V-6 4 1 



(5. 



-3/ 



Substituting (U) into (|5l|), we see that an expression of the form f = tih + t2h- X holds 
as long as one can find an X satisfying X -po = a*po, X -pi = api. Deriving A-po = apo, 
A ■ pi = a^pi from ( |5.3D and ( |5.7| ) , we obtain 

f = tih-t2h- A^ . (5.9) 

Notice that for fixed ti, the hi must be chosen such that the fi are integral. 



The relation (5.7) differs from that given in |5S] by the external factor a due to our different 



normalization for the periods ( |5.2| ). 
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The value of A^Aux, treating the degrees of freedom hi and tj as independent, is 



N^J^ = -t2h-A^-E-h, (5.10) 

where E is defined below (|2.2|). We see that A^Aux is independent of ti, indicating that 
the magnitude of the D3-brane charge does not care about the value of the RR axion 
Co = Re (j); this is consistent with the absence of a CqC^ A Hs A F3 coupling in the IIB 
supergravity action. 

The fraction of vacua at the LG point if; = can be estimated for Mq as follows. 
We can use the height-based method of analysis discussed in §2.4 to estimate the number 
of solutions to ( |5.9|) at fixed (j). Given the degree Vq = 2 and the fact that there is one 
complex structure modulus, we have ?7 = 4, reflecting the four real equations (|5.9|), so that 
from ( |2.2C1| ) we predict on the order of vacua for any fixed (p. More precisely, using 
( p.30|) we expect that for a given (p at height H{(p), 

iVvacua(i^;^,V' = 0)~|^ . (5.11) 

The factor of here just expresses the fact that roughly 1/i^ of the possible values 

for the hi are compatible with integral values of the fi on the LHS of each of the four 
equations in ( |5.9|) . We can therefore estimate the total number of vacua by summing over 
(/) G JFg and using the fact that there are 0{H'^) values of (/> at a given H, so that 

L 1 

iVvacua(^; = 0) ~ L2 ^ — . (5.12) 

H=l 

Note that the convergence of the sum over H at large L seems to indicate that the freedom 
to choose the dilaton ti, t2 does not contribute to the growth of the number of vacua; the 
scaling is consistent with the choice of the four hi alone. As emphasized in §2.4, this 
estimate does not take into account non-generic features of the constraint equations which 
may allow more or fewer solutions than the generic argument would predict. Indeed, as 
we see in the next subsection, there are indeed 0{L) solutions associated with (p at height 
L associated with special vacua where W = 0. The overall scaling of predicted by the 
generic height argument is still correct, however it misses this special class of vacua which 
contribute at large H and which also total to contribute order vacua. 
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While we have argued that the number of vacua at the LG point goes as L^, we 
expect that the total number of vacua over all of moduli space scales as ( |2.16| ) L^^ = L^. 
Therefore, the ratio of LG vacua to all vacua is 

j^^ . A^vacua(i^;V^ = 0) _ J_ 

The same result can be also argued from consideration of the restriction ( |5.1U| ) which 
for allowed fixed dilaton value gives the following scaling of the number of vacua 

N.acu.m = 0) ~ — ^ , (5.14) 

(Im (p)^ 

with an overall constant dependent on the arithmetic properties of the numerator and 
denominator of this scaling is in accord with the expected dilaton distribution in the 
fundamental domain. 

5.2.1 W = vacua 

Next we show that vacua with W = exist at the LG point for Mg. This equation 
leads to Vq = 2 additional constraints, and since the LG solutions (|5.9| ) still have four 
independent flixxes we generically expect solutions. The equation W = becomes 

{f-<j)h)-po = 0, (5.15) 

which is solved by LG flixxes satisfying 

h = h-A^ /i = (-3/i3 + /i4,3/i3,/i3,/i4) • (5.16) 

Together with ( |5.9| ) this implies / = / • A'^ as well. Hence requiring W — leads to a 
codimension 2 subspace in the lattice of fluxes. 

We can estimate the total number of W = vacua at the LG point using ( p.34|) and 
( p.35|) . One finds the total number of W = solutions will scale like 

iV,,eua(i^; ^ = 0; = 0) ~ 5^ ^ . (5.17) 

H=l 

Note that the sum over heights increases the L scaling by one power in this case, because 
of the modified contribution of the heights owing to the effective dilaton-independence of 
W = DfjjW = 0, as described in §2.5. 
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Another way to present these vacua is as follows. Instead of specifiying the fluxes 
/i3, /i4 and the dilaton parameters ti, ^2, one may instead specify fluxes /a, /4, /i3, /i4. The 
dilaton is then determined as: 

(a-2)/4-3/3 



with the D3-brane charge 



iVflux - 2(/3/i4 - Ms) • (5.19) 
We notice immediately the parallels with the generic vacua on the rigid CY studied in §3: 



four flux parameters determining the dilaton, as in (|3.5| ), and an analogous expression to 
(^) for iVfl,,. 

The counting of the vacua is hence identical to the rigid CY case. Fixing SL(2, Z)^ 
by choosing 

/4 = 0, 0</l3</3, (5.20) 



one obtains the scaling 



iVvacua(i^;V' = 0;Vr = 0) ~ — L^ (5.21) 

4o 



in agreement with ( 5.17 ). 

We turn now to studying the discrete symmetries of the model. As we shall see, an 
R-symmetry may be associated with the W — Q vacua, "explaining" the vanishing of the 
superpotential in the vacuum. 

5.2.2 Discrete symmetries on Mq 

The LG point is a flxed point for Zg C F, so one may hope that this symmetry is 
preserved in the low-energy theory. Additionally, the Z2 and Z3 points on the dilaton 
moduli space are also potential sources of low-energy symmetry. It is easy to see that 
given the requirement (|5.6| ) that the dilaton sits in jFg, the Z2 point is inaccessible, while 
vacua can sit on the Z3 point (ti = 0,ti = 1). The accessibility of the Z3 point is due to 
the sixth root of unity arising in both parts of the moduli space. 

We shall flnd in this case that the full flxed point symmetry Zg is not generically 
a symmetry of the low-energy action for LG vacua. However, there are cases where Z2, 
Z3 and even Zg symmetries do appear; the latter two, however, also involve a nontrivial 
Z3 C SL{2, Z). We shall start with the smaller symmetry groups and work our way up. 
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Z2 symmetries 

For W = vacua, we have thanks to (|5.16|) , 



W{cP, ^) = (/ - <ph) ■ n(^) = (/ - <j,h) ■ A' ■ n(^) = (/ - <ph) ■ (a'Uia^ij)) , (5.22) 

where we have used ( ^77|) . Since = —1, this imphes the superpotential is odd in i/j for 
all values of the dilaton: 

W{(l),-i;) = -W{(l),ij). (5.23) 

Hence a Z2 R-symmetry arises for these vacua, and it is easy to see that this forces 
VF(0, = 0) = 0. This is the R-symmetry "responsible" for the vanishing of the vacuum 



superpotential - in this example, this supports the suggestion in that R-symmetries 
enforcing the W = condition may not be too uncommon in the space of flux vacua. 
Clearly since "0 —>■ a^i(j, we can think of this Z2 as a subset of the full Zg acting only on 
moduli; the full Zg is not in general a symmetry. 

One may wonder whether analogous LG fluxes exist where the superpotential is even, 
and indeed this is the fluxes satisfying 

h^-h-A^ /i = (-3/13 + 3/i4, /i3, /i3, /i4) , (5.24) 

again together with the LG conditions ( p.9| ) imply 

W{(f),-ij) = W{(j),t(j) . (5.25) 

Hence for the fluxes ( p.24| ) we also have Z2 C Zg preserved, but as a true symmetry rather 
than an R-symmetry. No restriction is made on the value of the superpotential in this case 
and indeed it does not vanish. One can show instead that the vacua obeying (|5.24| ) satisfy 

D^D^W = 0. (5.26) 

Flux vacua satisfying DW = are imaginary self-dual and contain only (2, 1) and (0, 3) 
fluxes; W = is the condition for the vanishing of the (0, 3) part. Similarly, (|5.26|) is the 
condition that the (2, 1) flux vanishes.!! 



It is known that "attractor points of rank 2" can always be shown to host pure (0, 3) flux 



vacua; the LG point on Me is such a point |51,53]. 
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The counting of the D^D^W = vacua is entirely analogous to the W = Q case. 
Following the analysis of the previous subsection, we can either use the method of heights 
as in ( |5.17D , or we can note that given the expression for A^Aux for these vacua, 



iVfl 



ux 



(5.27) 



we again obtain the counting 




(5.28) 



Another perspective on these results is as follows. The expansion ( |5.2| ) for the periods 
indicates that the superpotential can be written in the form 



W{4>,^) = {f-<ph) ■ [poi^o(^') +Pimii^') +pli^'Fs{ij^) +pS^"^4(V'')] , (5.29) 



where we have used Ps = Pi,P4, = Po and C2+6n = c^+Gn = for all n, and the Fi are some 
Ze-invariant functions. As we have discussed, W = D^pW = imply f-po = h-po = f-pQ = 
h ■ Pq = 0, leaving VF((^, i/j) odd in i/j. Similarly, the combination D^W = D^^D^W = 
implies f ■ pi = h ■ pi = f ■ Pi = h ■ p\ =0, and consequently VF((/>, '0) for these vacua is 
even. 

Z3 symmetries 

Let us now study the possibility of the Z3 on the dilaton moduli space being preserved 
in the low-energy action. In order for fluxes to permit ^ = a as a vacuum, we must restrict 
the LG solutions ( |5.9| ) to 



Z3 is generated by the transformation cj) ^ {(p — 1) / (j). Is this a symmetry of the theory 
for the fluxes (HD? 

To see how this might work, recall the preservation of Z2 in the rigid case §3.5.1. 
There the fluxes transformed such that f — (ph {l/(p){f — (ph) ■ S, where S was a matrix 
of which n was an eigenvector. Similarly, it is possible to show that for fluxes ( |5.3C1| ), one 
finds under cf) {(j) — 1) / cp: 



f = -h-A\ 



(5.30) 



(/ - ^h) 



1 



(/ - 0/i)^' ■ 



(5.31) 
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Given (5]7), we find that the superpotential transforms as 

Vr(</>,^)^-Vr(</>,^), (5.32) 

provided we make the combined transformation 

(/)^ ((^- !)/(/), i)-^a^^. (5.33) 

Since ( |5.32|) indeed has the form W W/ (ccf) + d) up to the phase, this leads to a Kahler 
transformation ( |2.13| ) and hence a symmetry. Notice also that (just like the rigid Z2 case) 
when (j) takes its vacuum value, here (/> = a, we have W{(f) = a) — > Wi^cj) = a). 

Hence we find that for any LG solutions at the Z3 fixed point there is indeed a 
preserved Z3, but it is the diagonal between Z3 C SL(2, Z) and Z3 C Zg: the preserved 
symmetry of the low-energy action is a combination of both dilaton and complex structure 
modulus transformations. Without the latter, the former does not lead to a symmetry. 

The counting of such vacua can proceed as follows. Since we require that the dilaton 
sit at a fixed height, the number of vacua is simply determined by the counting of ( |2.20| ) 
with H{(f) = a) = 1, with the result 

iV,,,,,(L; = 0, (/) = a) ~ . (5.34) 

To make this more precise, note that for fixed dilaton (f) = a, the / fluxes are uniquely 
determined from h flux by ( |5.9|) with ii = and t2 = 1, and the number of vacua are then 
determined by positive quadratic form ( [5.1C1| ) again with t2 = 0. This gives 

iVvacua(i^;V' = 0,(^ = a) ~ — L2 . (5.35) 

So we see that discrete symmetries are generic at the LG point, and are down by at most 
a power compared to all vacua in this model. The fact that an order one fraction of 
all LG vacua sit at (p = a can be understood from the preference of the these vacua ( |5.11| ) 
to sit at small heights; since cp = a has i7 = 1, it is much more common among all LG 
vacua than some other dilaton with a larger height. 

Zg symmetries 

Now we at last consider vacua where the full Zq ofip —>■ aip can appear as a low-energy 
symmetry. The last example indicates that we may well need to combine this with some 
SL{2, Z) transformation on (j). Under ip — > aip we have 

W -^a*{f -(ph)-A- U{tlj) . (5.36) 
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Without an 5'L(2, Z) transform on we would require / and h both to separately be left 
eigenvectors of A; however, the eigenvectors of A are necessarily complex, and hence this 
is not possible. 

We find that with an SL(2, Z) transform on (p and certain choices of fluxes, a Zg 
symmetry can arise. In addition to the transformation 

1-0 



we must choose the fluxes so that both ( |5.32D as well as either ( |5.16| ) or (|5.24|) (in addition 



to the LG condition (|5.9| )) is satisfied. That is, is a good symmetry precisely when both 
Z3 and one of the Z2S is present: one needs the dilaton at (/> = ct as well as either W = 
or D(f)D^W = vacua. 



For the case of ( p.l6| ), the W = vacua with the Z2 R-symmetry, we find 

Wi<P,ij)-^——Wi<P,ij), (5.38) 



while for the D^j^D^W = vacua with non-R Z2 ( p.24| ) we have 

a* 

W{<P,ij)^^—^W{<P,ij). (5.39) 

These again have the proper 5'L(2, Z) form W — > W/ {c(f) + d) for a Kahler transformation 
( p.l3|) , and hence are symmetries of the theory. Moreover, it is simple to check that 
( ^.38|) and (|5.39| ) square to the Z3 transform (|5.32| ) , while they cube to the appropriate Z2 



transformations (|5.23| ), (|5.25| ). Hence, when both the Z3 and a Z2 symmetry are present 



simultaneously, they are subsumed into a full Zg. 

We finally count these two types of Zg vacua. From ( |2.34| ), the scaling for = Zg 
vacua at the fixed height H — 1 is simply given by 

iVvacua(i^; i; = 0,(P = a-W = 0)r^L, (5.40) 

and one can convince oneself that since D^W = Df^D^W = also become effectively 
dilaton-indepenent equations, the result for the D^pD^W = Zg vacua is the same. More 



precisely, one can note that the positive definite quadratic form (|5.10| ) for = solutions 
associated to ( |5.16|) reduces to 



iVflux = 2{3hl + 3hsh4 + hi) , (5.41) 
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while for the Dfj,D^W = solutions associated with (|5.24D it reduces to 



iVflux = 2{hl + hsh4 + hi) . (5.42) 
Consequently, in both cases the scaling becomes 



iVvacua(i^; i; = 0,(t> = a;W = 0)= iV,,,,,(L; ^ = 0, = a; D^D^W = 0) ~ 



TT 

2^3^ 

(5.43) 

Hence although the Z2 symmetries as well as the Z3 symmetry were generic within the set 
of LG vacua, all scaling as L^, we find that the intersection between the Z3 and one of the 
Z2 sets, which possesses a Zg symmetry, are supressed by 1/L. 

Complex conjugation C 

Before ending our discussion of the symmetries of Mq, we mention also the complex 
conjugation transformation C; our discussion is general and actually applies to all the 
Mfc. Using ( |5.2| ) one can demonstrate that the relation ( ^.38] ) indeed holds with Uk = 
rrik ■ V ■ fn'^^ , with the matrices ruk given in and 



V 



^0 IN^ 

10 

10 

Vi 0/ 



(5.44) 



This means that C indeed exists as a transformation carrying vacua to other (inequivalent) 
vacua: note that Afflux is preserved. 

A low-energy symmetry acting on the moduli alone will be present for fluxes invariant 
under C, that is f = f ■ Uk and h = —h ■ Uk- These vacua are present in all four models. 
For Mg, the conditions become 

h = -\h. /2 = -2/4, /l2 = /i3 = 0. (5.45) 

Not all of these vacua sit at the LG point; those that do additionally satisfy 

h = t2(2/i4 - hi) , /4 = t2{hi - ^/i4) , (5.46) 

and sit on the imaginary dilaton axis, = i^t2/2. In the set of LG vacua, we may 
think of those with the complex conjugation symmetry as those that also satisfy /i2 = 
/i3 = 0, ti = — 12/2. The overall scaling will go like L, but we should expect logarithmic 
enhancement along the lines of §3.5.2. We leave the question of the distribution of non-LG 
vacua with complex conjugation to the future. 
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5.3. Landau- Ginzhurg vacua for Mg 



The root of unity for the /c = 6 case was particularly "nice" because its real part is 
rational, leading to the small degree of the extension. The other cases are not so simple, and 
all three have extensions of degree four. The /c = 5 and /c = 10 cases both can be spanned 
with the bases {1, v^, za/S + y^, i\fh^ 5 + -\/5}, or alternately powers of the roots of unity 
such as {1, a, a^, a'*}. In general one may allow the dilaton = to + tia-\-t2a'^ + t4a^ with 
rational t^. Attempting to solve ( ^.4]) , one indeed finds eight equations, and no nonzero 
solution for the fluxes. So as we mentioned earlier, 



Ms, Mio : Ar,^cua(i^,V' = 0) = . 



(5.47) 



For the k = 8 case, the root of unity is a = exp(7rz/4) = (1 + i)/V^ with = i, and 
hence a convenient basis for this degree-four extension is {1, z, -\/2, z-\/2}. One considers 
the dilaton ansatz (j) = ti + it2 + \/2ts -\-i^/2t4, and solving the F-flatness conditions ( |5T^ ) 



with the vectors 



(8) 

Po 



nis- 



( ' \ 


/ 


a 




1 




V —ia / 


\ 



2 V2 
-1 



3 - i + v^- zV2 
-1 + 75 + 75 / 



(8) 
P2 



ms- 







ia 




1 




\-aJ 


\ 



- + - 



3 + i-V2-iV^ 

-1 - 75 + 75 / 
(5.48) 

one flnds in general that the fluxes are forced to vanish. However, for the simplifying 
assumption ts = t4 = something special happens: the four equations in D^W = 
become redundant with those in D(j,W = 0. 

One can see this as follows: the ifj equation and the conjugate of the (f) equation are 
{f — (ph) ■p2 = and {f -cph) -p^ = 0. It is evident from ( |5.48[ ) that p2 and Pq are identified 



under the exchange — V^, or equivalently a ^ a^; this is a transformation of the 

Galois group of jFg. Hence the F-flatness equations will be carried into each other under 
the same identiflcation, provided (p is invariant under \/2 — "\/2, which is precisely the 
requirement ts = = we found. Consequently, when cp contains no \/2 factors, the 
dilaton and ip conditions each reduce to the same four equations. 
The LG solution to these equations is: 



(f) = ti + it2 : 

with ti, t2 rational, and just as in the k — 6 case, 

f = tih- t2h ■ 



(5.49) 



(5.50) 
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with A this time given by [|59 |: 



A 



/I -10 1 ^ 

10-1 

-2-2 1 2 

V-4 4 1 -3/ 



(5.51) 



and also satisfying ( |5.7|) , but with the eighth root of unity. The expression for A^Aux in 
terms of the hi and ti is again 

^fluL = -t2h-A^-J:-h. (5.52) 

Notice that again the D3 charge is independent of the RR axion. 
As in Mg, one finds here that 

iVvacua(V' = Q)-L''J2jp- (5-53) 

This is because the counting works exactly for this nongeneric V = 4 case as it would in 
a generic T> = 2 case with one complex structure modulus: the sum over dilaton heights 
yields a factor of in (|2.31|) due to the specialization = t4 = 0, while there are only 



four independent linear constraining the flux integers instead of the expected eight. 

5.3.1 W=0 vacua 

Although this coincidence of the DW = equations leads to LG vacua where they 
might not have been expected, there is no similar coincidence for = vacua. Imposing 
the W = Q equation adds = 4 new conditions on the flux integers. Unlike the case for 
D(j,W = and D^W = 0, these conditions are not redundant with those already imposed. 
As a result = vacua must satisfy eight conditions on the eight fluxes, and the only 
solution is at vanishing flux integers. Hence there are no nontrivial W = Q vacua in this 
model. 

5.3.2 Discrete symmetries 

Consider now discrete symmetries. In addition to the Zg from the LG point, these 
models can access the Z2 point on the dilaton moduli space for ti = 0, t2 = 1? so these are 
the fixed points at which to look for low-energy symmetries; the Z3 point on the dilaton 
moduli space is forbidden. 

The most obvious symmetry comes from a property of the periods: as remarked before, 
all odd terms in tj) vanish in for the Mg and Mio models. Thus a Z2 is trivially present 
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under '0 — — i/j with W invariant. Although this symmetry has no action on the fluxes, it is 
nonetheless an element of the modular group, and in fact since a'^ = — 1, it is Z2 C Zg C F; 
modular transformations that do act on the fluxes square to this Z2. Hence, although it 
is a transformation of the moduli alone, it should be considered a gauge transformation 
of the theory, what we called in §2.6 a symmetry G, not a global transformation H; it is 
actually an identiflcation of the moduli space. 

Let us look for true global symmetries. We consider vacua at the dilaton Z2 fixed 
point, which must satisfy: 

f = -h-A^. (5.54) 
In a manner analogous to that of the k = 6 case, we can show that under (p — 1/^, 

{f-^h)^^{f-<f>h)-A\ (5.55) 

Consequently we have the combined transformation: 

(/>^ -!/(/>, V^^zV^, ^ Vr((/>,^)^^Vr((/>,^), (5.56) 

which is a Kahler transformation. Hence we flnd that the Z4 C Zg combines with the Z2 
of the dilaton to create what is naively a preserved Z4 symmetry in the low-energy theory. 
However, its square is the gauged Z2 described above; as a global symmetry, consequently, 
it is just a Z2. 

These vacua can be counted analogously to the Mq Z3 vacua in the previous subsec- 
tion; one takes (|5.53| ) and excludes the sum over heights, instead setting H{(j) = i) — 1: 

iVvacua(i^; (j) = i,lP = 0)r^L^ . (5.57) 

To obtain the coefficient, one again notes that for fixed dilaton ti = 0, t2 = 1 the / fiuxes 
are uniquely determined from the h fiuxes by ( |5.54| ) and the number of vacua is determined 
by the positive quadratic form ( |5.52|) . This gives 



7^2 

iVvacua(</' = Z, ^ = 0) ~ —L^ . (5.58) 

As with the Z3 vacua in Mg, we find that these Z2 vacua are generic among the LG vacua, 
both going as L^, because they sit at the lowest height for the dilaton. 

It is interesting to observe that we do not find any fiuxes with the full Zg preserved. 
A transformation —>■ aip with cj) inert would require / and h to be eigenvectors of A, but 
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these are complex. The Z4 (|5.56| ) is the most one can get from an SL{2, Z) transformation 
by S, and other SL{2, Z) transformations are less useful. 

As mentioned in the discussion of Mg, there is a complex conjugation operation C 
defined for all the M^, and for Mg there is a symmetry acting on the moduli alone for 
fluxes f = f ■ Us and h — —h ■ 11%: 

fi = -^fs, /2 = -2/4, h2 = hs=0, (5.59) 

with vacua at the LG point also satisfying 

f3 = t2{2h^ - hi) , h = t2{hi - 3h^) , (5.60) 

sitting on the imaginary dilaton axis, (p = it2] the counting is identical to the case for Mq. 

5.4- L-scaling of number of vacua near the LG and conifold points in Mg 

In (|5.2D and (pT^) we have presented the periods in an expansion about the LG point 



'0 = 0. We can use these expansions to perform a perturbative calculation of the vacua 
near the LG point. Although this has not been our primary focus in this paper, we give 
numerical results for the L scaling of the number of vacua in the vicinity of both the 
LG and the conifold points (an expansion of the relevant periods around the conifold can 
be found in [^), finding the expected scaling A^flux(-^) ~ = L'^ ( |2.16|) . Note that 



a similar calculation around the LG and conifold points in suffered from incomplete 
gauge fixing, which we correct in the present paper. This changes the scaling, bringing it 
into accord with the expected L^^ . 

In performing these Monte Carlo calculations, we looked for vacua at a variety of values 
-^flux ^ L. At a given A^Aux < i^, as in the analogous computations for the torus which 
were described in Section §4.3.2, we looked for vacua with fluxes max(|/|, \h\) < A;Aflux, 
where we used k = 1. Unlike for the torus, there do seem to be a nonzero fraction of vacua 
with fluxes outside this bound. Empirical tests, however, indicate that more than half 
of all fluxes at fixed Ag^x satisfy the bound in most of the region we considered. Thus, 
the scaling of the number of vacua as a function of L should not be affected, though the 
multiplicative constant we found may be too small by a factor of order 1. We could have 
imposed a larger bound, using, say /c = 2, which would capture almost all relevant vacua. 
This, however, would have slowed the program by a significant factor. 
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In both cases in our Monte Carlo computations we implement the following dilaton 
SL(2, Z) gauge fixing: 

/2=0, 0</li</i, /l27^0. (5.61) 

At special values of the fiuxes, additional discrete gauge identifications arise, but these do 
not change the L scaling. 

In the vicinity of the LG point, our numerical results give the following expected 
number of vacua for li/'p < 0.5 and N^^^ < L: 

N^acuaiL) = 0.02 L\ (5.62) 

Note that the LG Zg monodromy should be fixed in this case. We display our results in 
the figure below. Fixing the conifold monodromy by requiring arg(l — E [— tt, tt) one 
finds from Monte Carlo data the following expected number of vacua in the vicinity of the 
conifold point for |1 — -01 < 0.01 and A^Aux < L: 

N,acn.{L) = 0.2L\ (5.63) 

Our results for the expansion about the conifold point are also shown below. 

2 ■ lo"^ - 
1.5-10' - 

i-io' - 

500000 - 
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Fig. 7: The numerical results for the expected number of vacua in vicinity of LG 
point with iVfiux < L for L e (1, 100) fit by the curve 0.02 L". 
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20 40 60 80 100 

Fig. 8: The numerical results for the expected number of vacua in vicinity of 
conifold point with A^aux < L for L G (1, 100) fit by the curve 0.2 L^. 

6. Discussion 

In this paper, we have explored a variety of simple models for flux vacua. In particular, 
we have addressed the question of when flux vacua preserve supersymmetry in the no-scale 
sense, that is, which vacua have = 0, as well as flnding classes of flux vacua for which 
the low-energy fleld theories have discrete symmetries. Because in general the approach 
of approximating the fluxes by continuous parameters is not suited to solving the kind of 
overconstrained problems we encounter in studying these two questions, we have developed 
some simple approaches to addressing these types of questions, using techniques motivated 
from number theory. 

More speciflcally, we considered the models of the rigid Calabi-Yau, the symmetric 
torus with a single modular parameter, and the one-parameter Calabi-Yau models associ- 
ated with hypersurfaces in weighted projective space. For the flux vacua we considered, 
the periods (as well as their derivatives) live in a fleld JF of flnite extension degree T) over 
the rational numbers Q. In general, the number of constraints on the fluxes is linear 
in P, so that the form of the fleld extension plays an important role in determining the 
nimiber of allowed flux vacua. We found an exact solution for generic flux vacua on the 
one-parameter torus, with the result that the moduli of all such flux vacua live in a fleld 
extension of degree at most T> = Q over the rationals. On the one-parameter Calabi-Yaus, 
we considered vacua at the LG point, where the periods take values in the cyclotomic flelds 
arising from extending the rationals by k^^ roots of unity, with degrees T> = 2 and T> = A. 

In the models we considered, we found that the number of vacua with either 1^ = 
or with discrete symmetries is generally suppressed by a small power of L below the total 
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number of vacua, although in some cases we found anomafous scahng laws involving various 
powers of log L. These anomalous scalings arise from situations where the continuous 
approximation breaks down because either the geometry of the allowed volume in flux 
space includes a divergent region, or because extra constraints force us to use number- 
theoretic methods of analysis. 

For most of the models we considered, there is a single complex structure modulus and 
the total number of vacua goes as L^. We found that al\W = Q vacua on the symmetric 
torus have periods associated with degree 2 extensions over the rationals, and that the 
number of these vacua goes as L^. Families of vacua on the torus with discrete symmetries 
also generally scale as L^. In the one-parameter Calabi-Yau hypersurfaces that yielded 
vacua at the LG point, the results were similar: the total number of vacua was L^, and in 
fact the number of vacua with = or with discrete symmetries scaled with the same 
power of L as the total number of LG vacua, namely . In all these cases enhanced 
symmetry vacua are not suppressed by more than a power of relative to all vacua. 

While the precise results we have derived here are insufficient to make clear predictions 
for the scaling of = vacua and vacua with discrete symmetries when the number of 
complex structure moduli is large, we can gain substantial insight from the form of the 
equations we have used to do the analysis. In particular, as we have discussed here, 
imposing the equation W = ^ imposes an additional set of V equations on the fluxes at 
a point in moduli space where the periods live in a fleld extension of finite degree T> over 
Q. Since the number of additional constraints here (relative to generic vacua at the same 
point in moduli space) is not proportional to 63, we expect that generically when the total 
number of vacua (at the given point) grows as L"^, the number of = vacua will be on 
the order of L"~^/^, and will not be massively suppressed. 

Given the significance in our analysis of the arithmetic structure of the periods, an 
important question seems to be: in a Calabi-Yau with large 63, of the total of order L^^ 
fiux vacua, what fraction have periods which take values in fields of finite extension degree 
T) over the rationals? If these vacua with arithmetic periods form a large fraction of the 
total set of fiux vacua, then as long as P <^ 63, we expect that a significant fraction of 
the total number of vacua will have W = Q. Unfortunately it seems that currently little 
is known about this question. Based on the appearance of an unexpectedly low extension 
degree D < 6 for generic flux vacua on the torus, however, we are tempted to conjecture 
that the number of vacua with finite T> may be quite large. In fact, at this point we are 
unaware of any flux vacua which can be proven to have transcendental periods. A related 
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outstanding question is how to properly sum over families of vacua associated with different 
field extensions; more advanced concepts in algebraic number theory may be useful here. 

If the above conjectures are correct, so that indeed there are situations where W = 
vacua are suppressed by a power of L which is independent of 63, then 1^ = vacua may 
be "not too rare," in a way in which the authors of have suggested may indicate that a 
realistically small cosmological constant is not incompatible with low-scale supersymmetry 
breaking. Needless to say, a long and rather tenuous chain of reasoning relates the math- 
ematical problem we are addressing here to the phenomenological question of the scale of 
supersymmetry breaking. Numerous issues need to be understood much better before this 
connection can be convincingly justified. In particular, we need a better understanding of 
how a measure on the space of vacua can be produced by cosmological dynamics, and how 
such a measure can be incorporated into statistical methods of analysis such as were used 
2H] and in this paper. 
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Regarding discrete symmetries, based on the examples we have considered it seems 
likely that there are two types of discrete symmetries of flux vacua associated with Calabi- 
Yau manifolds with large 63. Some symmetries will involve imposing on the order of 63 
conditions on the fluxes, and should lead to a highly suppressed multiplicity of vacua, while 
other discrete symmetries may only require a fixed number of constraints independent of 
63, and may be much more prevalent. We analyzed symmetries in this paper which arise 
from discrete subgroups of the modular group as well as symmetries related to complex 
conjugation on the moduli space. The prevalence of these symmetries is related to the 
codimension of fixed point loci of the relevant group on the moduli space. Since the vacuum 
sitting at the fixed point locus is a necessary, but not sufficient condition, however, this 
codimension provides only an upper bound on the frequency of these vacua; a more exact 
understanding of under precisely which circumstances they arise would be of great interest 
for model-building and cosmology. 

There is a strong similarity between some of the questions we have investigated here, 
particularly the study of W = vacua, and the study of "rank two attr actors" associated 
with points on moduli space admitting fiux vacua with Gs of type (0,3) [ pT| , |5^ . In contrast, 
W = vacua arise at points in moduli space where Gs is purely of type (2,1), while the 
generic vacua satisfying the DW = equations are a combination of these two types. The 
attractor conjectures in their most recent form relate "rank two attractors" to Calabi-Yau 
varieties of CM type (see, e.g., for an explanation of this and related notions). In 
problems with /i^'^ = 1 (as we have mostly considered here) there is roughly a one-to-one 
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correspondence between those vacua with 6*3 of type (0, 3) and those with of type 
(2, 1). At large 63, however, it should be much more common to find points of the latter 
type than the former. While the (0,3) condition constrains to lie in a one- dimensional 
(complex) vector space, the W = Q condition leaves open an h?'^ dimensional subspace 
of , imposing as discussed above a number of additional constraints independent of 63. 
Thus, we expect that although there are similar issues involved in studying these two types 
of vacua, vacua with W = Q are likely to be much more prevalent than rank two attractors. 
The study of models with more than one complex structure parameter seems necessary to 
make progress on this issue. 

In this paper we have touched on potentially fascinating connections between the 
study of = vacua, vacua with discrete symmetries and number theory. For instance, 
we saw in §4.3.3 that the class numbers h{D) play an important role in the counting 
of = vacua on the torus, as conjectured in [^. More generally, we saw starting 
in §2.4 and running through all of our examples that simple ideas from the theory of 



heights [0 may prove to be of general use in exploring the taxonomy of string vacua. The 
stabilization of complex structure moduli through fluxes leads to a remarkable interplay 
between continuous geometry and discrete mathematics. It seems possible that further 
investigations of this connection will lead to a new and important role for number theory as 
a mathematical tool for understanding string compactifications and their phenomenology. 
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Appendix A. Equations for exact solution of symmetric torus 

As discussed in §4.3, the equations of motion DW = for the symmetric torus can 
be combined into a pair of quintic equations for the real and imaginary components of the 
modular parameter t = x + iy 

(li{x)y'^ =q3{x), 

(A.l) 

where 

qo{x) = acP — aPc , 

qi{x) = bc° - aPd + 2(ac° - aPc)x , 

q^(x) = boc - ado + (36c - 6oc° - Sad + a^do)x + (-36c° + 3a^d)x^ + 2(a°c - ac°)a;^ , 

q4{x) = —bod + bdo + 2{—boc + ado)x + (—36c + boc^ + Sad — a^do)x'^ 
+ 2(6c« - a°d)x3 + (-a«c + ac«)x^ . 

(A.2) 

These equations can be combined into a single cubic equation satisfied by x 

asx^ + a2X^ + aix + cu^ = , (A. 3) 

where the coefficients are given by 



q;° = a%lc^ - ablc^c^ - b^bo{c^fd + 2a%boc^d^ - {a^fbod^ - 2aa%oc^do 

+ 2aHocc^do + b^{c^fdo - 2a%'^c^ddo + {a^fbd^do + a^a^cdl - a^c^dl ; 



(A.4) 



ai = Qa%boc^ - Qabboc^c^ - 2a%lc^c^ - 2b%oc{c^f + 2ablc{c^f - 6aa%oc^d 
+ QaHocc^d + Sa^bbocc^d - Aabbo{c^fd - Q{a^fbocd^ + ^aa%oc^d^ 

- Qaa%c^do + 2{a^fboc^do + QaHcc^do - Aa%'^cc^do + Qab^{c^fdo (A.5) 

- 2aHo{c^fdQ + Qa^a^cddo + A{a^fbcddo - Ga^c^ddo - Saa^bc^ddo 
+ 2a{ayd''do - 2a{aycdl + 2a2a°c°(i^ ; 
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+ b^boicy - ablicy - ld>aa%c^d - 6{a^fboC^d + ISa'^bcc^d + Qa%^cc^d 

+ lQaa%QCC^d + ?,ab^{c^fd - 10a\{c^fd - 2a%bo{c^fd + 9a'^a^cd^ 

- Sia^bcd^ - da^c^d^ - 6aa%c^d^ + {a^fboc'd^ + Sa{a^fd' (A.6) 

+ 10{aybc^do - 16aa%cc^do - 2{aP)%occ^dQ + QaHic^do - a%\cydo 

+ Aa'^a^c^ddo + 2{a^)VddQ - {a^fd'^do + 2aa%Q{c^fdQ - 4a{a^fcddo 

+ {aycdl-a{ayc'dl, 

as = -Aia^boc^ - 6a%'^c^c° + ^aa%oc^c^ + 6a6^c(c°)^ - 4a^6oc(c°)^ 
+ 2a°66oc(c°)2 + 2b'^{c''f - 2abbo{cy + Gia^bc^d - 2{a'^)\cc'^d 
- GaHic^d - Ga^b'^ic^d + 2aa°6o(c°)^c^ - Gaia^cd'^ + Qa^a^c^d^ (A.7) 
+ 6(a°)26c°d2 _ 2(a°)^d^ + 4a{oPf^d^ - Sa^a^cc^do - 2{aybcc^do 
+ 4a^{cydo + 2aa°6(c°)^do + 2{aPfcddQ - 2a{aPfc^ddQ . 
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